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ON  THE  MAINTENANCE  OF  THE  POLAR  TRONT  JET  STREAM 


by 


J.  D.  Mahlman 


Abs  tract 


A  calculation  of  the  mean  transverse  circulation  about  the  polar  front 
jet  stream  is  performed  by  using  a  diagnostic  balance,  ui-equation  model. 
The  results  show  a  thermally-direct  mean  transverse  circulation  about  this 
jet  stream  system . 

An  examination  of  the  kinetic  energy  balance  of  this  jet  stream  reveals 
that  the  direct  transverse  circulation  is  strong  enough  to  maintain  the  jet 
against  frictional  dissipation  but  not  enough  to  provide  lateral  export  of 
energy.  However,  significant  amounts  of  energy  are  transferred  upward 
into  the  lower  stratosphere. 

Further  considerations  are  employed  to  show  that  the  mean  transverse 
circulation  obtained  here  is  compatible  with  the  observed  distributions  of 
temperature  and  potential  vorticity  about  the  jet  core. 


T 


1-2 


1 .  Introduction 


Since  its  quantitative  documentation  over  cwenty-five  years  ago,  the 
polar  front  or  mid-latitude  jet  stream  has  been  the  subject  of  numerous 
intensive  investigations.  Much  of  the  early  interest  in  this  jet  stream 
system  was  due  to  its  obvious  connection  with  the  mid-latitude  cyclone 
waves  and  their  associated  weather  processes.  Because  these  seemed 
to  form  an  interacting  system,  many  investigators  became  interested  in 
the  dynamical  processes  responsible  for  producing  and  maintaining  this 
jet  stream. 

According  to  Rossby  (1947)  and  Staff  Members  (1947),  the  westerly 
jet  stream  in  mid-latitudes  is  formed  by  lateral  mixing  processes  with 
constant  absolute  angular  momentum  to  the  south  of  the  jet  core  and  con¬ 
stant  absolute  vorticity  to  the  north.  This  assumption  yields  a  lateral 
velocity  profile  which  is  qualitatively  reasonable.  It  was  pointed  out, 
however,  that  such  a  mechanism  implies  a  thermally  indirect  transverse 
circulation  about  the  jet  core  itself.  Although  such  a  mechanism  is 
attractive  because  the  systematic  ascent  of  cold  air  and  descent,  of  warm 
air  helps  to  maintain  the  strong  lateral  temperature  contrast  across  the 
jet  core,  it  leads  to  the  possibly  more  difficult  question,  as  to  how  the 
kinetic  energy  of  the  jet  stream  is  maintained. 

Making  use  of  this  scheme,  Palmen  (1951)  included  this  assumed  in¬ 
direct  transverse  circulation  about  the  polar  front  jet  stream  in  his  model 
of  the  mean  meridional  circulation.  In  a  study  by  Endlich  (1953),  fields 
of  adiabatic  vertical  velocities  were  computed  in  the  vicinity  of  a  polar 
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front  jet  stream.  Although  fields  of  rising  and  sinking  motion  were  found 
on  both  sides  of  the  jet  core,  the  statistics  indicated  a  slight  preference 
for  the  indirect  transverse  circulation.  However,  no  systematic  averaging 
was  performed  relative  to  the  jet  core  itself.  In  a  separate  study,  Riehl 
and  Teweles  (1953)  indicated  results  which  were  not  incompatible  with  the 
original  concept  of  such  an  indirect  transverse  circulation.  Also,  Newton 
and  Carson  (1953)  and  Newton  (1958)  indicated  that,  in  the  vicinity  of  a 
developing  jet  stream  and  an  associated  frontal  zone,  the  vertical  motion 
branch  of  the  transverse  circulation  about  the  jet  core  must  be  in  an 
indirect  sense.  However,  they  contended  that  the  total  transverse  circu¬ 
lation  should  remain  thermally  direct. 

An  alternative  hypothesis  was  advanced  by  Namias  and  Clapp  (1949) 
who  proposed  that  the  jet  stream  forms  as  a  result  of  a  "confluence"  of  a 
preexisting  temperature  field.  In  this  case,  the  compensating  vertical 
transverse  circulation  is  thermally  direct.  A  similar  concept  was  also  pro¬ 
posed  by  Nyberg  (1949,  1950,  1953). 

The  thermally-direct  transverse  circulation  also  resulted  from  theore¬ 
tical  studies  by  Van  Miegham  (1950)  and  Kuo  (1954).  This  is  suggested  as 
well  in  theoretical  studies  of  atmospheric  frontal  structure  by  Sawyer  (1956) 
and  by  Eliassen  (1959,  1962),  Similar  results  were  indicated  in  a  numerical 
simulation  of  frontogenesis  by  Williams  (1967). 

The  only  previous  quantitative  attempt  to  establish  the  transverse 
circulation  about  the  polar  front  jet  stream  was  that  of  Clapp  and  Winston 
(1951).  They  found  tentative  evidence  in  favor  of  a  direct  circulation  by 
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using  adiabatic  vertical  velocities,  averaged  from  the  surface  to  10,000 
feet.  Similar  contentions  were  also  presented  as  a  result  ol  qualitative 
studies  by  Hubert  (1953),  Raethjen  (1953),  Vuorela  (1953  ,  !95  7),  Reiter 
(1961  ,  1963a),  Briggs  and  Roach  (1963),  and  by  Endlich  and  McLean  (1965). 

A  number  of  authors  have  investigated  the  vertical  motions  in  the 
vicinity  of  specific  portions  of  jet  streams.  Murray  and  Daniels  (1953), 
for  example,  showed  that  the  transverse  circulation  is  direct  in  th>  entrance 
zone  and  indirect  in  the  exit  zone  of  a  ’  igh  speed  center  along  the  jet  core. 

In  a  study  of  upper- tropospheric  frontogenesis ,  Reea  and  Sanders  (1953) 
showed  intense  sinking  motions  in  the  "jet  stream  front"  on  the  cyclonic 
side  of  the  jet  core.  This  was  corroborated  in  subsequent  studies  by  Reed 
(1955),  Reed  and  Damelsen  (1SS9),  and  Danielsen  (1959a).  These  authors 
argued  further  that  this  frontal  zone  may  contain  subsiding  stratospheric 
air  rather  than  being  a  zone  of  mixing  between  two  tropospheric  air  masses. 
This  contention  was  corroborated  in  case  studies  by  Danielsen  (1959b, 
lS64a,  b,  I960),  Danielsen,  Bergman,  and  Paulson  (1962),  Staley  (1960, 
1962),  Reiter  (1963a,  b,  1964),  Reitei  and  Mahlman  (1964,  1965a,  b) , 

Reiter,  Glasser  and  Mahlman  (1969)  and  by  Mahlman  (19b4,  1965).  However, 
this  intense  sinking  of  stratospheric  air  north  of  the  jet  core  was  shown  by 
Danielsen  (1964a)  and  by  Mahlman  (1964,  1965)  to  be  strongly  associated 
with  cyclogenesis  in  the  upper  troposphere.  Furthermore,  the  relationship 
of  this  exchange  mechanism  to  the  cyclogenetic  process  was  verified 
statistically  by  Mahlman  (1966,  1969). 

Even  though  the  above  mechanism  is  of  undeniable  importance  in  the 
exchange  of  mass  between  the  stratosphere  and  troposphere  and  in  the 
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dynamics  of  cyclone  development,  it  occurs  only  during  certain  periods 
along  limited  sections  of  the  jet  core.  As  a  result,  such  studies  may  not 
give  any  real  insight  into  the  nature  of  the  net  transverse  circulation  about 
the  polar  front  jet  stream . 

Prcbably  the  most  valuable  contribution  to  the  problem  of  the  mainten¬ 
ance  of  jet  stream  systems  was  that  presented  by  Riehl  and  Fultz  (1957, 
1958).  From  an  analysis  of  a  steady  three-wave  case  in  a  rotating-"dishpan" 
experiment,  the  authors  were  able  to  demonstrate  that  the  mean  transverse 
circulation  about  the  "dishpan"  jet  stream  is  thermally  direct.  This  was 
the  case  even  though  the  mean  meridional  circulation  in  the  same  latitudes 
was  determined  to  operate  in  a  thermally  indirect  sense.  Thus,  the  analysis 
showed  that  the  mechanisms  responsible  for  the  maintenance  of  the  jet 
stream  may  be  fundamentally  different  from  those  acting  to  maintain  the 
zonal-mean  wind  at  the  same  latitudes.  Riehl  and  Fultz’s  results  on  the 
energetics  of  the  “dishpan”  jet  stream  will  be  compared  later  against  the 
conclusions  of  the  analysis  to  be  performed  here. 

A  study  by  Elsberry  (1968)  of  a  dishpan  five-wave  vacillation  case  also 
indicated  the  existence  of  a  thermally-direct  transverse  circulation  around 
the  "dishpan"  jet  stream.  In  this  case  the  direct  circulation  was  a  sig¬ 
nificant  residual  between  thermally  direct  and  indirect  transverse  circu¬ 
lations  operating  at  the  poleward  and  equatorward  branches  of  the  jet 
stream  maximum,  respectively.  However,  in  the  cases  studied  by  Riehl  and 
Fultz  (1957,  1958)  and  by  Elsberry  (1968),  the  lateral  wind  shears  in  the  jet 
stream  regions  were  significantly  less  relative  to  the  Coriolis  parameter 
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than  those  observed  in  the  atmosphere.  Thus,  it  is  not  clear  whether  or  not 
such  results  can  be  immediately  applied  to  a  s/stem  such  as  the  polar  front 
jet  stream. 

The  success  of  Riehl  and  Fultz  (1957,  1958)  encouraged  efforts  to  obtain 
analogous  measurements  for  atmospheric  jet  streams.  Krishnamurti  (1961a, 
b)  in  a  detailed  study  of  the  subtropical  jet  stream ,  showed  the  existence  of 
a  thermally  direct  transverse  circulation  about  the  jet  core.  Thus,  the  trans¬ 
verse  circulation  about  the  subtropical  jet  stream  operates  in  the  same  sense 
as  the  mean  meridional  (Hadley)  circulation  in  the  same  region. 

In  a  study  of  the  stratospheric  polar  night  vortex,  Mahlman  (1966,  1969) 
found  a  direct  transverse  circulation  operating  below  the  polar  night  Jet 
stream  just  before  a  major  breakdown  period.  This  was  the  case  even 
though  the  mean  meridional  circulation  in  the  same  region  was  shown  to 
operate  in  the  indirect  sense . 

The  above  studies  have  been  useful  in  developing  our  understanding  of 
the  maintenance  of  jet  stream  systems.  The  polar  front  Jet  stream,  however, 
has  not  yet  been  studied  in  any  quantitative  detail  as  to  its  actual  trans¬ 
verse  circulation  and  corresponding  energetics.  This  is  to  a  great  extent 
due  to  the  difficulty  in  defining  an  adequate  curvilinear  coordinate  system 
for  such  a  strongly  time-dependent  phenomenon.  Investigators  have  been 
hampered  by  the  lack  of  diagnostic  techniques  capable  of  resolving  the 
vertical  motion  field  with  sufficient  accuracy  to  provide  reasonably  consis¬ 
tent  answers.  The  intent  of  this  study  is  to  obtain  a  systematic  calculation 
of  the  net  transverse  circulation  about  the  polar  front  jet  stream.  Further, 
it  is  desired  to  investigate  the  implications  of  such  a  circulation  on  the 
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energetics  of  this  jet  and  its  surroundings.  This  will  be  accomplished 
through  use  of  a  diagnostic  balance,  adequation  model  developed  originally 

j 

by  Krishnamurti  (1966,  1968). 

2 .  Experimental  Approach 

When  attempting  to  provide  a  meaningful  calculation  of  the  net  trans¬ 
verse  circulation  occurring  in  the  vicinity  of  a  polar-front  jet  stream, 
it  is  highly  desirable  that  the  case  selected  for  computation  be  one  in  which 
the  time  variations  of  the  system  are  not  particularly  large.  This  is  not 
only  advantageous  from  the  point  of  view  of  establishing  a  realistic 
coordinate  system ,  but  also  because  the  result  is  likely  to  be  more  repre¬ 
sentative  of  the  actual  processes  usually  acting  to  maintain  such  a  system. 

Because  of  the  above  considerations,  the  case  selected  for  analysis  was 
an  area  centered  over  the  contiguous  United  States  for  the  time  period 
15  November  1200  GMT  to  17  November  1200  GMT,  1966.  This  period  was 
dominated  by  a  polar  front  jet  stream  of  moderate  intensity.  The  flow  at 
jet-stream  level  was  somewhat  anticyclonic  over  the  continent  bounded  by 
relatively  weak  cyclones  off  each  coast  (Fig.  1).  Although  the  Pacific 
Coast  cyclone  did  move  slowly  eastward,  bringing  some  precipitation  to  the 
California  region,  no  significant  cyclongenesis  or  large-scale  system 
movements  were  observed  during  the  five  upper-air  observation  times  con¬ 
tained  in  this  selected  time  period.  Thus,  this  particular  case  appears 
satisfactory  for  attempting  to  determine  the  basic  dynamical  processes  acting 
to  maintain  the  polar  front  jet  stream. 
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Geopotential  height  (m)  at  300  mb  on  16  November  1966, 

0000  GMT,  Heavy  dashed  arrow  shows  position  of  the 
polar  front  Jet  stream  and  dashed  box  shows  lateral 
boundary  of  computational  grid.  °C  and  °A  represent  normal 
distance  in  degrees  latitude  on  the  cyclonic  and  anti- 
cyclonic  sides  of  the  jet  axis,  respectively. 


As  mentioned  in  the  Introduction,  the  approach  used  here  is  to  employ 
the  diagnostic  balance,  w-equation  model  developed  originally  by  Krishna- 

murti  (1966,  1968)  in  which  the  equations  are  scaled  for  motions  of  character¬ 
istic  Rossby  number  less  than  unity  (R  <  1).  The  initial  data  input  is 
obtained  from  the  National  Meteorological  Center  objectively  analyzed  geo¬ 
potential  fields  at  1000,  850,  700,  500,  300,  and  200  mb.  The  stream 
function  (^)  is  assumed  to  be  related  to  the  geopotential  field  through  the 
non-linear  balance  equation,  * 

7.fV0+  2J  ,  |^)=  V2*  (1) 

4> 

In  this  study  £  is  specified  at  the  lateral  boundaries  by  simple  ^  . 

o 

In  the  event  that  the  eilipticity  condition  for  Eq.  (1), 

v24»+jf2-7f*V|/)>0,  (2) 

is  not  satisfied,  Eq.  (1)  is  replaced  by 

v  •  fv^  =  V2#  -  2J  (ug,  vg)  •  (3) 

Such  an  adjustment  is  frequently  required  on  the  anticyclonic  side  of  the 
jet  stream.  Use  of  Eq.  (2)  in  such  cases  in  place  of  Eq.  (1)  permits 
retention  of  strong  anticyclongic  vorticities  from  the  diagnosed  height  fields . 

Ey  using  the  balance  equation  (I)  combined  with  the  continuity,  vor- 
ticity,  and  thermodynamic  equations,  one  obtains  the  "complete" 

(Ro  <  1)  «-equation  system, 

^See  the  Appendix  for  the  list  of  symbols. 
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repeated  until  the  difference  between  successive  approximations  of  the  w 
field  become  acceptably  small.  In  this  study,  five  passes  were  found  to  be 
sufficient.  For  further  details  in  the  computational  proc  dure,  see  Krlsh- 
namurti  (1968). 

The  computations  are  made  from  a  five-level  grid  in  the  vertical. 

Fields  of  0  ,  % ,  4>  ,  8^/dt  appear  at  the  1000,  800,  600,  400,  and  200  mb 

surfaces  while  w,  T,  9  ,  q,  and  q  are  given  at  the  900,  700,  500,  and 

s 

300  mb  levels.  The  horizontal  grid  points  are  spaced  2.5  degrees  latitude 
and  longitude  apart.  The  grid  extends  from  55W  to  135W  longitude  and  from 
25 N  to  60N  latitude.  There  are  IS  grid  points  in  the  meridional  direction 
and  33  grid  points  in  the  zonal  direction.  The  last  six  grid  points  do  not 
contain  initial  data,  but  are  used  to  provide  cyclic  continuity  so  that  grid 
points  1  and  33  have  the  same  value  for  any  given  dependent  variable.  Grid 
point  values  for  points  26  through  32  are  provided  by  interpolation  from 
values  at  grid  points  1,  2,  26,  and  27. 

Once  the  calculations  for  each  observation  time  are  completed,  a 
coordinate  line  is  defined  along  the  polar  front  Jet  stream  axis  at  300 
mb  (as  defined  from  a  careful  isotach  analysis).  Then  parallel  coordinate 
lines  are  defined  at  distances  2,5  and  5.0  degrees  latitude  perpendicular 
to  the  jet  axis  in  the  cyclonic  and  anticyclonic  directions,  respectively 
(see  Fig.  2).  It.  is  then  a  comparatively  straightforward  procedure  to 
calculate  "mean"  distributions  of  variables  relative  to  the  jet  core.  For 
an  example  of  the  mean  wind  and  temperature  structure  relative  to  this  jet, 
see  Fig.  2.  In  this  system,  an  "eddy"  quantity  is  defined  as  the  value  of 
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a  variable  at  a  given  point  minus  its  "mean"  along  the  same  line  parallel 
the  jet  axis.  At  times  it  will  be  useful  to  discuss  "mean"  and  "eddy" 
processes  relative  to  the  jet  core.  These  are  to  be  carefully  distinguished 
from  the  mean  and  eddy  concepts  arising  from  the  zonal  averaging  often 
employed  in  general  circulation  studies. 

3 .  Mean  Transverse  Circulation  About  the  Jet  Core 

The  calculations  of  uj,  averaged  along  one  complete  wavelength  of  the 
polar  front  Jet  stream  were  prepared  for  all  levels  at  the  specified  normal 
distances  to  the  300  mb  jet  core.  These  were  prepared  individually  for 
each  observation  time  from  15  November  1200  GMT  to  17  November  1200 
GMT,  1966. 

The  individual  calculations,  derived  from  the  five  observation  times, 
all  showed  for  the  troposphere  a  distinct  thermally  direct  pattern  of  mean 
vertical  motion  about  the  jet  core  (i.e. ,  systematic  ascent  of  the  warm 
air  on  the  anticyclonic  side  of  the  jet  core  relative  to  the  colder  air  on  the 
cyclonic  side).  Although  each  independent  calculation  revealed  this,  in 
the  interest  of  increasing  the  reliability  of  subsequent  inferences,  the 
individual  calculations  were  combined  to  form  a  time  averaged  mean  uc 
relative  to  the  jet  core. 

The  quantitative  results  of  this  are  given  in  Fig.  3  along  with  the  95 
percent  confidence  interval,  determined  from  the  time  variability  of  the 
individual  measurements.  This  figure  shows  a  systematic  mean  ascent  on 
the  anticyclonic  side  of  the  jet  core  with  maximum  values  at  300  mb.  This 
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Fig.  3.  Mean  vertical  motion  ([u))s>  10  mb  sec  )  relative 

to  the  jet  axis  averaged  over  the  period  15  November 
1200  GMT  to  17  November  1200  GMT,  1966.  Numbers  to 
the  right  of  the  points  give  the  95^  confidence 
interval  determined  from  the  calculations. 


characteristic  agrees  well  with  observations  of  cirrus  cloud  shields  on  the 
anticyclonic  side  of  the  jet  core  (McLean,  1957;  Conover,  1960;  Kadlec, 

1963;  Oliver,  Anderson,  and  Ferguson,  1964;  Bittner,  1966;  Anderson, 
et  al, ,  1969).  On  the  other  hand.  Fig.  3  shows  ascending  motion  through 
the  jet  core.  This  result  may  disagree  with  satellite  studies  which  contend 
that  the  cirrus  shield  cuts-off  at  the  jet  core  or  just  to  the  anticyclonic  side 
(Oliver,  Anderson,  and  Ferguson,  1964;  Bittner,  1966;  Anderson,  et  al. , 

1969).  Nevertheless,  this  calculated  ascent  in  the  jet  core  is  not  incon¬ 
sistent  with  the  result  of  Starrett  (1949),  which  shows  that  precipitation  tends 
to  be  a  maximum  directly  underneath  the  jet  core.  It  is  also  compatible 
with  the  distributions  of  ozone  and  water  vapor  mixing  ratio  in  the  vicinity 
of  the  jet  core,  given  by  Briggs  and  Roach  (1963).  This  result  is  further 
supported  by  the  observation  that  the  tropopause  "gap"  (apparently  dynam¬ 
ically  produced)  is  generally  located  just  to  the  cyclonic  side  of  the  jet 
core . 

In  a  study  such  as  this,  where  a  number  of  types  of  averaging  procedures 

must  necessarily  appear,  it  is  convenient  to  employ  the  averaging  notation 

proposed  by  Reiter  (1969a,  b).  In  this  system  the  mean  of  an  arbitrary 

quantity,  say  Q,  is  given  by  [Q ] ^ ,  where  the  quantity  enclosed  by  small 

parentheses  indicates  the  coordinate  (or  coordinates)  over  which  the 

averaging  was  performed.  The  deviation  of  Q  from  that  average  is  given 

by  (Q),..  Thus,  for  example,  one  might  write  Q  =  [Qj,  .  +  (Q).  .  for 
0  (s)  (s) 

averaging  processes  along  the  jet  axis. 


A  useful  way  to  view  the  mean  motion  relative  to  the  jet  core  is 
through  use  of  a  mass  transport  stream  function.  Because  of  the  large  ratio 
of  the  length  to  the  width  of  the  Jet  stream  In  this  study,  one  may  write  the 
averaged  mass  continuity  equation  as 


w  _  0 

an  ap  0  ' 


(7) 


One  can  write  a  solution  to  this  equation  in  the  form  of  a  stream  function 
such  that 


[vj]e> 


ara,.,  »w,., 

- £si  tvjj  =  +-—1M 

ap  '  ^(s  an 


(8) 


The  diagnostic  balance,  co-equation  model  is  solved  with  a  boundary  condi¬ 
tion  of  w  =  0  at  100  mb.  Because  of  this,  ^  is  arbitrarily  set  equal  to 
±ero  on  that  boundary.  The  implied  net  inflow  into  the  region  is  obtained  by 
integrating  Eq.  (7)  from  5  degrees  latitude  on  the  anticyclonic  side  to  5 
degrees  latitude  to  the  cyclonic  side  of  the  jet  core.  It  is  then  assumed  that 
the  net  computed  inflow  or  outflow  can  be  partitioned  equally  on  the  anti¬ 
cyclonic  and  the  cyclonic  boundaries.  Once  the  lateral  boundary  values 
for  [Vj]^sj  are  set,  then  the  entire  field  of  ^  is  determined  by  numerical 
integration  of  Eq.  (8).  The  results  of  this  determination  along  with  the  mean 
wind  structure  of  the  region  are  given  in  Fig.  4.  This  figure  shows  a  number 
of  interesting  features.  Rather  strong  confluence  of  warm  and  cold  air  appears 
to  be  taking  place  at  about  the  700  mb  level  under  the  jet  core.  This  is  in 
good  agreement  with  the  concept  p-esented  by  Namias  and  Clapp  (1949). 

Note  the  well-defined  flow  across  the  jet  axis  at  200  mb  and  the  horizontal 
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splitting  of  streamlines  on  the  anticyclonic  side  of  the  jet  axis.  This  implies 

rather  strong  energy  generation  through  -  [v_]  .  .  ~  [♦],  .  at  the  jet  core. 

j  (s )  a  n  ts; 

Also,  there  is  an  implied  horizontal  divergence  on  the  anticyclonic  side  of 
the  jet  core.  This  is  compatible  with  the  observed  lateral  spreading  of 
cloud  elements  in  the  cirrus  shield  on  the  anticyclonic  side  of  the  subtropical 
jet  stream  as  seen  in  ATS-I  satellite  time-lapse  movies.  Further,  a  study 
by  Fujita,  et  al.  (1968)  of  these  ATS-I  time  lapse  movies  shows  significant 
lateral  convergence  of  plumes  from  the  tops  of  lower-  to  mid-tropospheric 
convective  cells  to  the  south  of  the  cirrus  shield.  This  is  also  in  agreement 
with  the  results  given  in  Fig.  4. 

Ic  is  of  interest  to  note  that  the  structure  given  in  Fig.  4  is  also  quite 
similar  in  overall  character  to  the  qualitative  model  of  transverse  circulation 
about  the  polar  front  jet  core  given  by  Danielsen  (1968)  for  cases  of  strong 
import  of  stratospheric  air  into  the  "jet  stream  front"  in  association  with 
cyclogenesis.  This  appears  to  be  of  fundamental  interest,  since  the  case 
selected  for  analysis  here  had  no  occurrences  of  large-scale  cyclogenesis 
during  the  period.  Therefore,  even  though  the  sinking  of  stratospheric  air 
on  the  cyclonic  side  of  the  jet  core  in  such  cases  is  about  an  order  of 
magnitude  more  intense  than  the  mean  values  given  here  (see  numerous 
references  on  this  phenomenon  cited  in  the  Introduction),  the  qualitative 
forms  of  the  transverse  circulations  are  very  similar.  Because  of  this 
correspondence,  it  is  then  reasonable  to  hypothesize  that  the  fundamental 
dynamical  mechanisms  responsible  for  maintenance  of  the  polar  front  j€.-t 
stream  and  for  producing  rapid  upper  tropospheric  cyclogenesis  (with 
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associated  import  of  stratospheric  air  into  the  troposphere},  may  be  essentially 
the  same.  This  may  be  the  case  even  though  the  former  exhibits  a  quasi¬ 
steady-state  character  while  the  latter  is  essentially  discrete  in  nature. 

This  hypothesis,  of  course,  needs  much  more  thorough  examination  before  its 
validity  car.  be  strongly  claimed.  Such  comparisons  and  examinations  are 
difficult,  not  only  because  these  are  energetically  open  systems,  but  also 
because  they  interact  very  strongly  with  each  other. 


4.  Jet  Stream  "Mean"  Kinetic  Energy  Maintenance 

Given  the  mean  transverse  circulations  about  the  polar  front  Jet  stream 
from  the  previous  section,  it  is  possible  to  provide  quantitative  estimates 
of  the  various  contributions  which  this  circulation  produces  toward  maintaining 
the  "mean"  kinetic  energy  of  this  jet  stream.  The  conventional  approach  to 
evaluate  the  energetics  of  open  systems  for  volumes  which  do  not  encircle 
the  globe,  is  to  derive  expressions  for  the  time  rate  of  change  of  the  total 
kinetic  energy  of  the  given  volume.  Since  the  major  Interest  here  is  in  the 
maintenance  of  the  "mean 11  structure  of  the  polar  front  jet  stream,  only  the 
balance  equation  tor  the  mean  kinetic  energy  (averaged  along  the  jet  axis) 
will  be  evaluated  in  any  detail. 

The  derivation  of  the  mean  kinetic  energy  balance  in  an  arbitrary 
curvilinear  (s,n,p)  coordinate  system  is  obtained  by  first  considering  the 
quasi-static  equations  of  motion  in  this  system, 


6  u  _  3u  3u  3u  o$ 

—— 1  +  u  — -i-  +  v  — -L  -f  oj— ■ -  fv  +  —  +  F  =  0  , 
3 1  J  3s  J3n  3p  J3s  s 


(9) 
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*VT  9VI  3VS  3VI  34 

*r~  +  u  rJ  v  — —  +  w  -rJ-  +  fur  +  —  +  F  =  0  .  (10) 

3t  J  3s  J  3n  3p  J  3n  n 

It  should  be  pointed  out  here  that  the  s  axis  Is  along  the  jet  axis  and  not 
necessarily  along  the  instantaneous  streamlines.  Also,  is  the  wind 
component  along  this  jet  axis  and  v  Is  the  wind  component  along  the  normal 
n  axis . 

Now,  by  multiplying  £q.  (9)  by  [u  and  Eq,  (10)  by  [Vj]^,  adding , 
averaging  along  s  and  then  along  n,  manipulating,  and  integrating  over  the 
mass  contained  between  two  pressure  surfaces,  one  obtains  the  equation  for 
the  time  rate  of  change  of  "mean"  kinetic  energy  of  the  Jet  stream. 


->  <- 


+ 


+  J  A  „  Cr 


(n)  +  Vw'Vb)  +  [vI](s)(vJ>fe)  dt 
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In  the  above  expression  the  separate  terms,  labeled  (a)  through  (j) 
have  the  following  interpretation: 

(a)  flux  of  "meaif  kinetic  energy  through  the  side  boundaries; 

(b) ,(c)  flux  of  "mean"  kinetic  energy  through  the  lower  and 
upper  boundaries,  respectively; 

(d)  conversion  of  "eddy"  kinetic  energy  to  "mean"  kinetic 
energy; 

(e)  an  additional  term  in  the  conversion  from  "eddy"  to  "mean" 
kinetic  energy  which  arises  because  the  Coriolis  parameter 

—  (f)  varies  along  the  axis  over  which  the  initial  average  is 
performed.  It  is  interpreted  as  a  conversion  term  because 
it  appears  with  opposite  algebraic  sign  in  the  "eddy" 
kinetic  energy  equation; 

(f)  mean  pressure  interaction  term  at  the  side  boundaries  arising 
from  energy  conversion  in  an  open  system; 

(g) ,(h)  mean  pressure  interaction  terms  at  the  lower  and  upper 
boundaries,  respectively; 

(i)  conversion  of  potential  energy  into  "mean "  kinetic  energy  due 
to  mean  transverse  circulation  about  the  jet  core; 

0)  dissipation  of  "mean"  kinetic  energy  due  to  friction  and  sub¬ 
grid-scale  mixing. 

Because  the  coordinate  system  tends  to  be  aligned  nearly  parallel  to 
the  flow,  not  only  are  the  terms  involving  "eddy"  products  difficult  to 
calculate,  but  many  times  they  may  prove  to  be  negligibly  small.  Since  the 
intent  of  this  work  is  to  assess  the  effect  of  the  "mean"  transverse 
circulation  about  the  Jet  stream,  only  quantities  aJf  b^  c  ,  f,  g,  h,  and 
i  are  computed  here.  A  program  to  provide  explicit  calculation  of  "eddy" 
effects  is  being  presently  developed.  As  will  be  pointed  out  later,  there 
m  ly  be  some  aspects  of  jet  stream  dynamics  which  are  strongly  dependent 
upon  "eddy"  processes  relative  to  the  jet  stream.  For  the  present  analysis, 
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however,  the  "mean"  kinetic  energy  balance  is  put  in  the  simplified  sym¬ 
bolism  implied  in  Eq.  (11), 

— -  =  a  +  b,  +  c  +f+g+h+i  -  "dissipation  .  "  (12) 

o  t  I  1  1 

Here,  "dis  Uution"  is  actually  a  combination  of  terms  b^,  c 2<  and 

j.  As  a  result,  its  algebraic  sign,  as  well  as  its  magnitude,  is  presently 

uncertain.  The  approach  here  is  to  assess  the  values  for  the  remaining 

terms  as  well  as  for  3K.,/dt  so  that  some  knowledge  can  be  gained  about  the 

M 

combined  "eddy"  and  dissipation  effects.  The  calculations  are  performed 
from  900  to  100  mb  over  an  area  five  degrees  latitude  perpendicular  to 
the  jet  core.  In  order  to  gain  further  clarification  of  the  processes  involved, 
separate  calculations  are  performed  for  the  volumes  100  -  300  mb  and  300  - 
900  mb.  This  roughly  approximates  the  regions  above  and  below  the  jet 
core. 

The  results  of  these  calculations  are  given  .n  Table  1.  First,  noting 

the  results  for  the  entire  100  -  900  mb  volume,  it  may  be  seen  that  terms 

a  ^ ,  b^,  and  in  Eqs.  (11)  and  (12)  produce  a  small  net  loss  of  energy 
-2  -1 

(-0.8  erg  cm  mb  sec  )  to  the  volume  due  essentially  to  mean  lateral 

outflow  (see  Fig.  5)  occurring  at  krger  mean  kinetic  energies  than  the 

mean  inflow.  The  mean  pressure  interaction  terms  (f,  g,  and  h  in  Eqs.  (11) 

and  (12)  also  show  a  net  loss  for  the  100  -  900  mb  volume  during  this  time 
-2  -1  -1 

period  (-1.8  erg  cm  mb  sec  ).  The  conversion  from  potential  to  mean 
kinetic  energy  of  the  jet  (term  (i)  in  Eqs.  (11)  and  (12))  shows  a  net  gain 
of  +1.9  units  for  the  entire  volume.  This  conversion  is  especially  significant 
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Table  1 
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since  the  local  change  term  (SK^/dt)  represents  a  net  loss  of  5.1  units  for 


the  period.  The  net  transverse  circulation  about  the  polar  front  Jet  is  thus 

shown  to  be  thermally  direct  even  when  the  jet  is  rapidly  becoming  less 

intense.  As  a  result,  the  mean  transverse  circulation  calculated  here  may 

be  less  intense  than  in  cases  where  the  jet  may  be  in  approximate  steady- 

state  or  accelerating.  The  implied  "dissipation"  required  for  balance  in 

-2  -1  -1 

this  case  is  4.6  erg  cm  mb  sec 

A  number  of  studies  have  been  undertaken  recently  by  Kung  (1966a,  b; 

1967;  1968;  1969)  on  the  generation  and  dissipation  of  kinetic  energy  in  the 

atmosphere.  It  is  interesting  to  note  that  his  most  recent  estimate  of  the 

-2  -1 

annual  average  energy  dissipation  over  North  America  is  4.7  erg  cm  mb 
sec  1  for  the  100  -  900  mb  layer.  Although  this  value  is  virtually  identical 
with  the  implied  "dissipation"  obtained  here,  there  are  a  number  of  important 
differences.  The  transfer  of  kinetic  energy,  due  to  the  unevaluated  eddy 
terms  in  Eq.  (11),  is  potentially  important  and  could  be  of  either  algebraic 
sign.  Also,  Kung's  dissipation  values  are  for  the  total  kinetic  energy  and 
not  for  the  mean  jet  stream  component,  as  calculated  here.  It  is  usually 
assumed  that  the  energy  dissipation  of  the  largest  scales  of  motion  is 
appreciably  smaller  than  for  smaller-scale  motions  (e.g.,  see  Oort  (1964)). 

An  alternative  hypothesis  has  been  advanced  by  Panofsky,  (1970)  who 
speculated  that  since  the  vertical  scale  of  ruesoscale  and  larger-scale  systems 
is  not  appreciably  different  (particularly  so  in  the  k '  stream  region),  large- 
scale  kinetic  energy  might  be  transferred  directly  into  sub-grid-scale  motions 
without  going  through  any  significant  scale  cascading  processes.  Some 
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evidence  for  this  is  given  by  Clarke  (1966)  who  estimated  a  mean  dissipation 
-2  -1  -1 

rate  of  4.5  erg  cm  mb  sec  for  the  100  -  500  mb  layer  in  the  vicinity 
of  the  subtropical  jet  stream.  If  this  is  a  valid  statement,  it  suggests  that 
a  great  deal  of  the  estimated  "dissipation"  in  this  study  may  be  true 
dissipation.  This  suggests  that  only  a  comparatively  small  amount  of  the 
kinetic  energy  produced  is  exported  out  of  the  jet  stream  region.  IYesently, 
this  can  be  only  speculation  at  best,  because  of  combined  effects  of 
computational  uncertainty,  unevaluated  "eddy"  terms,  and  the  observed 
deceleration  of  this  Jet.  It  is,  however,  in  agreement  with  the  results  of 
Rrehl  and  Fultz  (1958),  who  stated  for  the  three-wave  steady  "dishpan" 

Jet:  "One  cannot  view  the  jet  stream  as  a  current  that  is  the  locus  of 
generation  of  kinetic  energy,  which  is  then  dispersed  laterally  to  maintain 
the  kinetic  energy  elsewhere,  nor  as  a  current  that  is  'driven  by  outside 
circulations , . . . " 

If  one  partitions  the  calculations  into  balances  for  the  regions  below 

(300  -  900  mb)  and  above  (100  -  300  mb)  the  jet,  considerable  additional 

information  is  gained.  For  the  region  below  the  jet,  Table  1  shows  a  very 

pronounced  decrease  of  mean  kinetic  energy  with  time  ,/bt  =  -10.6  erg 

M 

—2  -i  -1 

cm  mb  A  sec  .  It  is  especially  interesting  that  the  conversion  to  mean 
kinetic  energy  (i)  is  +3.0  units.  There  are  significant  losses  through  the 
boundary  terms  (c1  and  h)  across  the  300-mb  surface.  For  this  volume  the 
implied  "dissipation"  is  6.3  units.  As  before,  it  is  not  presently  possible  to 
separate  out  the  eddy  terms  from  the  true  energy  dissipation.  Given  reason¬ 
able  estimates  of  the  dissipation  (Kung's  (1969)  value  for  the  900  -  300  mb 
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layer  is  5.8  erg  cm  mb  sec  )  and  the  validity  of  Panofsky's  (1970) 
hypothesis,  the  unevaluated  eddy  terms  in  Eq.  (11)  may  not  be  particularly 
important  for  this  volume.  This,  however,  remains  to  be  tested 
quantitatively. 

The  region  above  the  jet  appears  to  possess  entirely  different  character¬ 
istics.  Table  1  shows  that  this  layer  is  increasing  its  kinetic  energy 

-2  -1  -1 

(3KM/at  =  +11.2  erg  cm  mb  sec  ).  A  large  portion  of  this  increase 
is  due  to  flux  of  mean  kinetic  energy  (term  b)  into  the  region  across  the 
300  mb  surface  (+9.6  units).  Also,  a  very  large  flux  of  energy  into  the  region 
is  produced  by  the  mean  pressure  interaction  term  at  300  mb  (+14.3  units). 
Some  of  this  energy  is  lost  through  the  side  boundaries  (term  a^81  -3.9  units 
and  term  f  =  -8.0  units).  It  is  especially  interesting  to  note  that  the  conver¬ 
sion  term  (i)  is  a  -1.3  units  and  thus  acting  to  build  up  potential  energy  at 
the  expense  of  mean  kinetic  energy.  This  arises  because  of  the  reversal  of 
the  mean  lateral  temperature  contrast  above  the  jet  core.  This  is  also 
compatible  with  the  contention  of  Riehl  (1962)  that  this  temperature  reversal 
above  the  jet  core  is  produced  dynamically  by  ascent  on  the  anticyclonic  side 
and  descent  on  the  cyclonic  side  of  the  jet  core.  The  implied  "dissipation" 
for  this  layer  is  -0.5  units.  This  value  is  too  small  to  be  considered  reliable, 
even  with  respect  to  algebraic  sign. 

A  truly  steady-state  condition  would  have  been  more  favorable  for  inter¬ 
pretative  purposes.  It  is  quite  encouraging,  however,  that  the  summation  of 
various  kinetic  energy  balance  terms  are  quite  compatible  with  the  observed 
tendencies.  This  suggests  that  the  diagnostic  balance,  oj-equation  model 
is  useful  for  such  studies.  It  does  seem  clear  that  the  mean  polar  front  jet 
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stream  is  maintained  by  a  thermally  direct  transverse  circulation  below  the 
core,  and  that  appreciable  energy  is  transported  upward  into  the  lower 
stratosphere  indirectly  as  a  result  of  this  process. 

These  results  carry  some  interesting  implications.  First,  it  is  clear 

that  the  mechanisms  for  maintaining  the  jet  stream  at  these  latitudes  are 

completely  different  from  the  processes  acting  to  maintain  the  zonal  mean 

westerlies,  as  viewed  in  the  traditional  approach  to  general  circulation 

problems.  Therefore,  it  appears  to  be  inconsistent  to  equate  these  two 

phenomena  in  any  way  when  discussing  the  mechanisms  involved  in  their 

maintenance.  The  second  implication  is  that  the  polar  front  jet  stream 

provides  a  very  efficient  mechanism  for  producing  upward  flux  of  energy  into 

the  stratosphere  due  to  the  pressure  interaction  effect.  In  fact,  the  flux 

across  the  300  mb  surface  due  to  term  (g)  in  Eq.  (11)  and  Table  1  is  2860 
-2  -1 

erg  cm  sec  .  In  the  usual  general  circulation  formulation  of  the  energy 
equations,  this  would  show  up  as  a  very  strong  eddy  effect  at  about  wave 
numbers  4  or  5 .  The  above  value  is  quite  compatible  with  the  vertical  eddy 
pressure  interaction  term  for  the  same  general  region  obtained  by  Manabe  and 
Hunt  (1968)  in  an  18-level  general  circulation  model.  It  is  especially  interest 
ing  to  note  that  this  large  energy  flux  into  the  stratosphere  occurs  at  the 
approximate  location  of  the  winter  mid-latitude  warm  belt  in  the  lower  strato¬ 
sphere.  This  region  is  known  to  be  a  strong  sink  of  zonal  and  eddy  available 
potential  energy.  Thus,  the  above  process  may  be  of  fundamental  importance 
in  maintaining  the  kinetic  energy  of  this  somewhat  unusual  region  of  the 
straosphere.  It  would  be  of  interest  to  check  this  hypothesis  directly  through 


a  high-resolution  numerical  model  of  the  general  circulation. 

5.  Some  Additional  Considerations 

In  the  previous  section  considerable  attention  was  given  to  the  balance 
of  "mean"  kinetic  energy  along  the  polar  front  Jet  stream.  One  can  also 
write  a  companion  expression  to  Eq .  11  for  the  "eddy"  kinetic  energy 
balance  along  the  jet  axis. 

In  the  present  experimental  setup,  it  has  been  impossible  to  calculate 
all  terms  in  such  an  "eddy"  kinetic  energy  balance  because  of  the  difficulty 
in  obtaining  "eddy"  products.  However,  as  a  means  of  obtaining  a  prelim¬ 
inary  estimate  of  the  energetic  contribution  of  "eddies"  relative  to  the  polar 
front  jet  stream,  a  calculation  of  the 

term  in  the  "eddy"  kinetic  energy  balance  was  performed.  Physically,  this 

represents  a  conversion  from  potential  energy  into  "eddy"  kinetic  energy 

measured  relative  to  the  jet  axis.  This  particular  term  is  quite  important 

because  it  helps  to  determine  whether  the  "eddies,"  relative  to  the  jet  axis, 

have  an  internal  energy  source  or  have  been  produced  by  either  boundary 

fluxes  or  conversions  from  "mean"  kinetic  energy.  Figure  5  shows  the 

average  of  the  above  expression  over  the  five  observation  times  at  each 

grid  point  (in  terms  of  energy  units).  The  volume  average  of  this  conversion 
-2  -1  -1 

is  +5.1  erg  cm  mb  sec  .  This  represents  a  rather  strong  energy  source 
in  the  jet  region  and  is  probably  large  enough  to  balance  the  sub-grid  scale 
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Distribution  of  potential  into  "eddy"  kinetic 
energy  conversion  relative  to  the  Jet  axis 
averaged  over  the  period  15  November  1200  GMT 
to  17  November  1200  GMT,  1966.  Units  are 
expressed  in  erg  cm  ‘•mb  sec  . 


energy  dissipation  of  "eddy"  kinetic  energy  without  necessitating  a  large 
input  from  the  various  other  possible  mechanisms. 

Figure  £  shows  that  this  energy  conversion  has  an  interesting  distri¬ 
bution  relative  to  the  jet  axis .  A  very  strong  source  of  "eddy"  kinetic 

-2  -1  -1 

energy  (~20  erg  cm  mb  sec  )  exists  in  the  middle  troposphere  on  the 
cyclonic  side  of  the  Jet  axis.  This  source  is  apparently  due  to  the  asymmetric 
distribution  of  the  cyclone-scale  waves  relative  to  the  jet  axis.  On  the 
anticyclonic  side  of  the  jet  axis  in  the  lower  troposphere,  a  significant 
(but  weak)  sink  of  energy  appears .  This  strong  variation  across  the  jet 
axis  is  somewhat  surprising.  Its  dynamical  significance  is  not  presently 
obvious,  although  it  seemingly  must  play  an  important  role  in  the  redistri¬ 
bution  of  kinetic  energy  in  the  jet  stream  region.  The  effect  appears  to  be 
statistically  significant  because  it  appears  in  each  of  the  five  calculations 
at  the  individual  observation  times.  As  in  any  case  study,  however, 
companion  calculations  for  other  cases  would  yield  valuable  insight  into 
the  generality  of  these  preliminary  results  . 

As  pointed  out  by  Riehl  (1962)  and  others,  one  of  the  principle  diffi¬ 
culties  with  a  thermally  direct,  transverse  circulation  about  the  jet  axis  is 
that  such  a  circulation  acts  to  rotate  the  isentropic  surfaces  into  horizontal 
planes.  Thus,  the  solenoidal  field  necessary  to  support  a  geostrophic  jet 
stream  is  weakened  by  the  transverse  circulation.  Yet,  in  the  first  approxi¬ 
mation,  geostrophic  balance  holds  for  the  Jet  stream.  This  Implies  a 
restriction  on  the  intensity  of  the  transverse  circulation  depending  upon  the 
intensity  of  the  processes  capable  of  counteracting  the  adiabatic  cooling  on 
the  anticyclonic  side  and  the  adiabatic  heating  on  the  cyclonic  side  implied 
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in  Fig.  3.  Some  insight  can  be  gained  into  the  possible  processes  respon¬ 
sible  by  considering  the  formulation  for  the  mean  heat  balance  along  the 
jet  axis.  By  averaging  the  first  law  of  thermodynamics  along  the  jet  axis 
and  separating  "mean"  and  "eddy"  processes  one  obtains 


(«)  (b) 


A  preliminary  investigation  has  been  made  of  the  possible  processes 
represented  in  Eq.  (12)  which  may  counteract  the  anticyclonic  side  cooling 
and  cyclonic  side  heating  resulting  from  term  (f)  in  the  same  equation.  This 
analysis  indicates  that  the  only  term  which  can  accomplish  this  is  (b),  which 
represents  the  convergence  of  "eddy"  heat  flux  relative  to  the  jet  axis.  In 
order  for  this  to  be  valid,  an  "eddy"  heat  flux  divergence  is  required  on  the 
cyclonic  side  and  an  "eddy"  convergence  on  the  anticyclonic  side.  This 
conclusion  was  also  reached  by  Riehl  and  Fultz  (1958)  in  their  analysis  of  a 
three-wave  "dishpan"  jet.  As  indicated  previously,  a  program  is  now  under¬ 
way  to  provide  for  explicit  calculation  of  such  "eddy"  terms. 
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Another  very  important  feature  of  the  structure  of  an  atmospheric  jet 

stream  is  that  the  potential  vorticity  (P  =  -  ■—  (f  +  £  J)  is  almost  discon- 

oP  o 

tinuous  across  the  jet  axis.  As  pointed  out  by  Riehl  (1962)  such  a  potential 
vorticity  distribution  in  the  vicinity  of  the,  jet  axis  places  a  strong  restriction 
on  the  type  of  dynamical  process  that  can  be  responsible  for  maintenance 
of  the  jet  stream.  Realizing  this,  It  may  be  instructive  to  consider  whether 
or  not  the  transverse  circulation,  alleged  here  to  be  responsible  for  main¬ 
tenance  of  the  polar  front  jet  stream,  is  compatible  with  the  observed  strong 
gradients  of  potential  vorticity  across  the  Jet  axis.  In  this  case  a  potential 
vorticity  equation  considering  non-conservative  effects  is  appropriate,  since 
the  time  required  for  a  parcel  to  make  a  complete  circuit  about  the  Jet  axis 
(as  obtained  from  Fig.  4)  is  the  order  of  40-50  days.  As  shown  by  Staley 
(1960),  potential  vorticity  can  only  be  modified  by  a  curl  of  the  frictional 
force  and  gradients  of  diabatic  heating.  (It  should  be  pointed  out,  however, 
that  such  a  parcel  would  not  in  general  make  such  a  complete  circuit  about 
the  axis  because  of  rather  large  "eddy"  and  sub-grid  scale  turbulence 
effects .) 

A  preliminary  analysis  was  performed  to  estimate  possible  magnitudes  of 
the  various  terms  in  the  "mean"  potential  vorticity  balance  equation.  The 
transverse  circulation  of  Fig .  4  and  reasonable  estimates  on  the  friction  and 
diabatic  heating  were  employed  to  provide  these  estimates.  The  analysis 
showed  that  these  non-conservative  terms  are  important  on  the  longer  time 
scales  but  are  not  directly  responsible  for  producing  the  observed  potential 
vorticity  gradient  across  the  jet  axis.  The  results  show  that  the  probable 
cause  for  this  is  that  the  jet  core  occurs  at  the  tropopause  level,  and  the  net 


1-34 


transverse  circulation,  as  seen  in  Figs.  3  and  4,  leads  naturally  to  such  a 
potential  vorticity  distribution.  On  the  cyclonic  side  of  the  jet  axis,  stable 
stratospheric  air  is  being  slowly  brought  downward  while  tropospheric  air 
is  being  slowly  lifted  on  the  anticyclonic  side.  Because  of  the  relatively 
small  net  cross-axis  velocity  component  required  to  maintain  the  jet  ener¬ 
getically,  the  mean  advection  normal  to  the  jet  core  tends  to  be  considerably 

smaller.  This  is  particularly  true  at  the  jet  core  itself.  In  more  quantitative 

3[P](s) 

terms ,  the  -  [(»>],  .  — - — *-*■  term  in  a  potential  vorticity  balance  relative 

(S)  P  ^(S) 

to  the  jet  axis  clearly  dominates  the  -  [vj]^  — — 1  term.  Although  the 
horizontal  convergence  of  "eddy"  potential  vorticity  flux  cannot  yet  be 
evaluated  directly,  it  is  probable  that  this  effect  is  not  strong  across  the 
jet  core  itself,  because  streamlines  at  this  level  tend  to  be  quite  parallel 
to  the  axis  of  maximum  wind. 

Furthermore,  from  the  case  studies  of  discrete  intrusions  of  stratospheric 
air  into  the  troposphere  mentioned  in  the  Introduction,  considerable  vertical 
"eddy"  flux  of  potential  vorticity  is  to  be  expected  in  the  vicinity  of  the 
polar  front  jet  stream.  In  this  case  the  resultant  vertical  convergence  of 
"eddy"  potential  vorticity  flux  will  also  act  to  reinforce  the  observed  near 
discontinuity  across  the  jet  axis.  The  essential  similiarity  of  this  process 
(which  occurs  only  at  particular  regions  and  times  along  the  jet  axis)  to  the 
net  transverse  circulation  has  already  been  pointed  out  in  section  3.  As  a 
result,  the  conclusion,  that  these  are  reinforcing  effects,  is  not  surprising. 

Thus,  one  may  conclude  that  the  deduced  net  transverse  circulation  about 
the  jet  core  is  not  contradicted  by  cither  heat  balance  requirements  or 


potential  vorticity  considerations 
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6.  Summary  and  Conclusions 

Calculations  have  been  made  of  the  mean  transverse  circulation  about 
the  polar  front  jet  stream,  using  a  diagnostic  balance,  w -equation  model. 
The  results  show  that  such  a  circulation  is  thermally  direct  in  the  tropo¬ 
sphere,  with  ascent  on  the  anticyclonic  side  and  through  the  jet  core  and 
descent  on  the  cyclonic  side. 

A  calculation  of  the  "mean"  kinetic  energy  of  the  jet  stream  shows  that 
the  polar  front  jet  is  probably  maintained  by  this  direct  circulation,  and 
that  no  appreciable  amounts  of  energy  are  either  imported  into  or  exported 
out  of  the  jet  stream  region.  The  results  do  show,  however,  the  existence 
of  a  comparatively  large  upward  flux  of  energy  into  the  stratosphere, 
resulting  from  the  thermally-direct  transverse  circulation. 

Preliminary  analyses  of  the  heat  and  the  potential  vorticity  budgets 
indicate  that  the  mean  transverse  circulation  obtained  here  is  compatible 
with  the  observed  temperature  and  potential  vorticity  distributions  in  the 
vicinity  of  the  polar  front  jet  stream. 
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APPENDIX 


Table  of  Symbols 

A  area  enclosed  by  lateral  boundaries 

B  index  denoting  averaging  or  integration  to  be  performed  along 

lateral  boundaries 

c  velocity  component  normal  to  the  lateral  boundary  (+  inward) 

n 

Cp  specific  heat  of  air  at  constant  pressure 

f  Coriolis  parameter 

fg  average  f  over  region 

Fs<Fn  frictional  force  in  the  s  and  n  directions,  respectively 

g  acceleration  of  gravity 

H  diabatic  heating  rate 

Hl  diabatic  heating  rate  due  to  release  of  latent  heat 
J  Jacobian  operator 

Km  mean  kinetic  energy  measured  with  respect  to  the  jet  stream  axis 
1  length  along  the  boundary  B 

n  coordinate  direction  oriented  perpendicular  to  the  jet  stream  axis 

P1,Pu  pressure  at  lower  and  upper  boundaries  of  integration  volume, 
respectively 

q,qg  specific  humidity  and  saturation  specific  humidity,  respectively 
R  gas  constant  for  dry  air 

Rq  Rossby  number 

s  coordinate  direction  oriented  along  the  jet  stream  axis 

T  temperature 

Uy,Vg  zonal  and  meridional  components  of  geostrophic  wind,  respectively 
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urvf  Wind  components  along  and  normal  to  the  Jet  stream  axis,  respectively 
Vg  wind  speed 

ft  specific  volume 

8  Rossby  parameter 

C'C  8  relative  vorticity  on  p  and  9  surfaces,  respectively 
C3  t  +  r 
x  R/c 

P 

7  "  dp  '  a  static  stability  parameter 

Tx'Ty  frictional  stress  on  the  x  and  y  directions ,  respectively 
$  geopotential 

X  velocity  potential 

^  stream  function  for  horizontal  wind 

*  aST™  fUnCti°n  f°r  mean  frans verse  circulation  about  Jet  steam 

tc  dp/dt 
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ABSTRACT 


A  numerical  model  which  employs  observed  stratospheric  winds  to 
advect  simulated  tracers  was  developed.  This  model  successfully  repro¬ 
duces  many  qualitative  features  of  the  observed  fields  of  both  ozone  and 
radioactive  debris.  As  the  tracers  evolve,  the  horizontal  eddies  con¬ 
stitute  the  principal  process  modifying  the  tracer  zonal  mean.  Although 
the  vertical  eddies  and  the  zonal  mean  cell  are  both  an  order  of  magnitude 
weaker  than  this  process,  the  latter  of  these  tends  to  always  act  in  the 
same  sense  so  that  its  effect  becomes  more  important  over  long  periods 
of  time  . 
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I.  BACKGROUND 


A.  OBSERVATIONS 

The  study  of  trace  substance  transport  Ln  the  stratosphere  can  pro¬ 
vide  many  useful  insights  into  the  dynamics  of  that  region.  The  history 
of  an  evolving  inert  tracer  depends  upon  its  initial  state  and  upon  several 
statistical  properties  of  the  wind  field.  In  the  stratosphere,  unprocessed 
observations  of  the  wind  field  are  not  always  adequate  to  reveal  these 
rather  subtle  characteristics  of  the  general  circulation.  In  addition  to 
observations  of  the  wind,  there  exist  some  observations  of  the  distri¬ 
butions  of  such  natural  tracers  as  ozone  and  radioactive  debris.  These 
observations  may  be  used  to  extend  the  available  data  base,  and  theories 
of  the  general  circulation  may  be  tested  by  requiring  them  to  simulate 
the  advection  of  natural  tracers . 

The  influence  of  stratospheric  transport  processes  was  first  noticed 
when  the  early  spectrophotometer  observations  of  Dobson  and  his  co¬ 
workers  (Dobson  et  al.  1928)  revealed  that  the  distribution  of  total  atmos¬ 
pheric  ozone  has  the  following  characteristics: 

1 .  The  total  amount  of  ozone  per  unit  area  increases  from  the  equator 
toward  the  poles. 

2.  In  high  latitudes  this  quantity  varies  with  season,  exhibiting  a 
pronounced  maximum  in  late  winter  or  early  spring. 

3.  No  apparent  correlation  exists  between  total  ozone  and  solar 
activity . 
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4.  Total  ozone  is  correlated  with  the  passage  of  upper-level  synoptic 
systems,  attaining  positive  departures  from  its  long-term  mean  in  the 
lows  and  negative  departures  in  the  highs. 

Dobson's  observations  are  consistent  with  those  obtained  by  the  81 
station  network  in  the  Soviet  Union  and  reported  by  Bozkov  (1968).  These 
observations  Indicate  that: 

1.  Minimum  total  ozone  occurs  at  10  N. 

2.  Maximum  total  ozone  occurs  at  60  to  70  N,  with  values  decreasing 
toward  both  the  equator  and  the  pole. 

3.  The  maximum  zonal  mean  gradient  of  total  ozone  occurs  between 
30  and  60  N. 

4.  The  total  ozone  tends  to  have  a  wavelike  distribution  in  the 
horizontal  with  low  values  occurring  over  the  western  portions  of  conti¬ 
nents  and  high  values  over  the  eastern  portions. 

The  photochemical  theory  of  ozone  formation  (See  Craig,  1965)  pre¬ 
dicts  that,  at  photochemical  equilibrium,  the  distribution  of  total  ozone 
should  attain  a  maximum  in  low  latitudes  at  the  subpolar  point  and 
increase  toward  the  poles.  Brewer  (1949)  and  Dobson  (1956)  attempt  to 
reconcile  this  theory  with  observation  by  postulating  a  mean  cell  with 
ascending  motion  in  low  latitudes,  northward  advection  in  the  high  mid¬ 
latitude  stratosphere,  descent  near  the  pole,  and  a  return  circulation  at 
low  levels  in  mid-latitudes.  This  circulation  supposedly  carries  the 
ozone,  which  is  generated  by  solar  ultraviolet  radiation  in  low  latitudes, 
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northward  and  downward  into  the  lower  polar  stratosphere  where  it  ac¬ 
cumulates  since  there  it  is  protected  from  photodissociation. 

■j 

The  Brewer-Dobson  theory  is  seriously  challenged  by  airborne  filter 
observation  of  the  distribution  of  Tungsten-185  from  the  Hardtack  nuclear 
tests  (Feely  and  Spar,  1960).  The  center  of  concentration  of  this  debris 
remained  at  50  mb  and  nearly  at  the  latitude  of  injection.  The  cloud  did 
not  migrate  northward  and  downward,  but  rather  spread  in  the  meridional 
plane  in  such  a  manner  that  the  isopleths  of  concentration  sloped  down¬ 
ward  from  the  equator  toward  the  pole . 

B.  PARAMETERIZED  MODELS 

Reed  (1950)  and  Normand  (1953)  note  that  stratospheric  troughs  are 
regions  of  both  high  temperature  and  high  total  ozone.  This  leads  them 
to  postulate  that,  since  both  potential  temperature  and  ozone  mixing 
ratio  increase  upward,  the  lows  must  be  regions  of  descending  motion 
and  the  highs  regions  of  ascent.  Furthermore,  the  troughs  arri  ridges 
m<-)uo  mona  cim»/jv  the  v.rl nc!  so  thst  horlzcr.ts!  ^dvGcticr.  zYi if t c  th»c 
centers  of  greatest  mixing-ratio  departure  downwind  from  the  associated 
pressure  system.  Because  the  wind  almost  always  has  a  westerly  com¬ 
ponent  in  the  winter  season,  the  ozone  maxima  lie  in  the  southerly  winds 
cast  cf  the  troughs  and  the  minima  are  east  of  the  ridges  in  northerly 
winds.  Thus,  die  departures  of  both  temperature  and  mixing-ratio  are 
explained  m  terms  of  the  systematic  interaction  between  horizontal  and 
vertical  a  r.ection.  High  mixing-ratio  and  temperature  are  accompanied 
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by  descent  and  poleward  motion  while  low  values  of  these  two  quantities 
occur  with  rising  and  equatorward  motion. 

Citing  earlier  work  by  Dobson  (1960),  Newell  (1961)  carries  this 
line  of  reasoning  one  step  further.  He  explains  the  poleward  slope  of 
the  isopleths  by  invoking  eddy  mixing  in  which  northward  motion  is 
correlated  with  descent  nd  southward  motion  with  ascent.  In  the  same 
paper  he  computes  the  flux  of  ozone  by  correlating  total  ozone  with  the 
winds  at  12-18  km  in  the  maximum  ozone  layer.  He  finds  that  90  percent 
of  his  computed  flux  is  due  to  large-scale  eddies  and  that  the  total  flux 
is  more  than  adequate  to  explain  observed  changes  in  the  distribution. 

Prabbakara  (1961)  generates  zonal  nean  ozone  distributions  which 
depend  upon  the  interaction  of  photochemistry,  a  mean  cell  that  descends 
near  the  pole,  and  large-scale  eddies  parameterized  as  anisotropic  dif¬ 
fusion.  His  results,  which  are  in  good  qualitative  agreement  with  obser¬ 
vation,  show  that  such  transport  processes  can  explain  many  departures 

C - —  -I - i  —  -  1  ---.**.•»  ■ 

v»i.»  (#.I  1UC  A  1U1II  * 

In  contrast  to  Prabhakara's  work,  in  which  the  principal  axis  of 
eddy  diffusion  is  horizontal,  there  have  been  a  number  of  models  simu¬ 
lating  the  transport  of  radioactive  debris  with  anisotropic  diffusion  schemes 
in  which  the  principal  axis  is  inclined  to  the  horizontal  (Reed  and  German, 
1965;  Davidson,  et  al.  ,1965;  Seitz,  et  al.  ,1968;  and  Fairhall  and  Peed, 
1966).  In  such  models  it  is  assumed  that  the  northward-downward  velo¬ 
city  correlation  in  the  large-scale  eddies  can  be  parameterized  in  terms  of 
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a  diffusion  tensor  whose  principal  axis  slopes  downward  toward  the 
pole.  This  means  that  the  tracer  is  presented  with  a  path  of  least  re¬ 
sistance  parallel  to  the  observed  slope  of  the  isopleths  of  mixing  ratio 
in  the  real  atmosphere.  In  order  to  keep  the  tracer  on  these  sloping  paths 
and  to  prevent  excessive  vertical  spreading,  it  is  necessary  to  make  the 
vertical  component  of  the  diffusion  tensor  much  smaller  than  that  along 

the  principal  axis.  Davidson  et  al.,  (1966)  used  a  principal  diffusivlty 

9  2  -1  4  2  -1 

of  4  x  10  cm  sec  and  a  vertical  diffusivity  of  about  10  cm  sec 

However,  since  characteristic  horizontal  distances  are  on  the  order  of 

3 

10  km  while  vertical  distances  are  about  1  km  and  the  speed  of 
diffusion  depends  upon  the  square  of  these  distances,  a  ratio  of  106 
between  the  two  diffuslvities  is  not  physically  unreasonable. 

C.  GENERAL  CIRCULATION  MODELS 

Hunt  and  Manabe  (1968)  and  Hunt  (1969)  report  the  most  extensive 
and  realistic  simulation  of  stratospheric  tracers  to  date.  In  their  experi¬ 
ments  various  simulated  substances  are  introduced  into  the  wind  fields 
generated  by  a  general  circulation  model,  and  the  history  of  the  numeri¬ 
cally  evolved  mixing  ratio  fields  is  examined. 

In  the  1968  experiments  inert  tracers  are  permitted  to  evolve  from 
zonally  symmetric  initial  states.  The  first  such  tracer  roughly  corre¬ 
sponds  to  radioactive  debris  from  an  equatorial  injection  and  is  initially 
distributed  in  a  band  extending  from  the  equator  to  ten  degrees  north 
with  the  maximum  concentration  at  50  mb  pressure  height. 
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The  second  experiment's  initial  state  represents  the  distribution 


of  ozone  at  photochemical  equilibrium  throughout  the  region  of  integration. 
Both  tracers  are  treated  as  inert  substances  inasmuch  that,  except  for 
removal  of  any  tracer  that  finds  its  way  into  the  troposphere,  no  sources 
or  sinks  are  included  in  these  experiments. 

The  simulated  histories  of  both  tracers  are  reassuringly  similar  to 
observations  of  real  tracers  in  the  atmosphere.  Both  the  radioactive 
dobris  and  the  ozone  begin  at  once  to  migrate  northward  and,  to  some 
extent,  downward.  Thus,  the  concentration  at  the  equator  decreases  and, 
at  the  poles  increases,  so  that  the  initial  gradient  becomes  smaller  and 
eventually  reverses. 

The  actual  transfer  by  the  quasi-horizontal  eddies  is  well  illustrated 
by  the  evolution  of  the  mixing  ratio  field  at  various  horizontal  levels  in 
the  model.  The  zonally  symmetric  character  of  the  initial  field  quickly 
disappears ,  becoming  sinusoidal  with  wave  number  four  predominating  as 
it  does  in  the  circulation  used  by  the  model.  With  the  passage  of 
time  the  waves  becomes  increasingly  exaggerated,  forming  the  mixing- 
ratio  contours  into  long  northward  protrusions.  Eventually  the  protrusions 
Separate,  leaving  isolated  islands  of  high  concentration  in  northern  lati¬ 
tudes  . 

Hunt  (1969)  reports  two  additional  experiments  with  the  same  model 
and  involving  ozone  transport  with  photochemical  sources  and  sinks  added. 
The  mechanism  of  transfer  in  these  studies  and  in  those  previously 
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discussed  are  very  similar.  In  both,  the  mean  cell,  which  produces  a 
convergence  of  tracer  into  the  subtropics,  is  opposed  by  an  o  Idy  diver¬ 
gence  out  of  that  region,  with  the  eddies  being  particularly  effective 
in  transferring  tracer  northward  in  the  region  poleward  of  30N.  The  prin¬ 
cipal  difference,  introduced  by  photochemistry,  is  the  strong  source  in 
low  latitudes  which  prevents  reduction  or  reversal  of  the  initial  poleward 
gradient  by  replacing  the  ozone  as  fast  as  transport  processes  can  re¬ 
move  it. 

The  weak  link  in  any  experiment,  involving  either  parameterization 
or  a  general  circulation  model,  lies  in  the  possibly  unrealistic  specification 
of  the  transport  processes,  and  its,  at  best,  indirect  relation  to  the  actual 
motions  in  the  stratosphere.  In  view  of  these  difficulties,  it  seems  instruc¬ 
tive  to  compare  advection  of  simulated  tracers  by  actual  observed  winds 
with  both  observation  and  previous  attempts  to  model  both  ozone  and  radio¬ 
active  debris.  The  present  study,  therefore,  employs  actual  observed 
stratospheric  winds  to  advect  simulated  tracers  in  the  presence  of  sub¬ 
grid-scale  diffusion,  and  numerically  integrate?  the  equation  of  continuity 
for  the  tracer  in  time  to  obtain  the  evolution  of  an  inert  tracer  from  an 
arbitrary  initial  state. 

This  approach  is  largely  independent  of  any  preconceived  ideas  about 
the  stratospheric  general  circulation,  and  so  provides  a  third  source  of 
information  intermediate  between  the  parameterized  models  or  the  general 
circulation  models  and  observations  of  real  tracers. 
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II.  MODEL 


A.  BASIC  EQUATION 

The  present  study  follows  the  lead  of  Hunt  and  Manabe  (1968)  by 
numerically  exploring  the  detailed  advection  of  trace  substances  without 
resort  to  arbitrarily  defined  diffusion  tensors.  Furthermore,  the  winds 
used  to  accomplish  this  advection  are  based  on  actual  observations  of 
the  real  atmosphere  and  are  not  the  result  of  a  mathematical  model  that 
may  or  may  not  reproduce  the  properties  of  the  actual  general  circulation. 

The  basic  equation,  on  which  this  study  depends,  is  the  equation  of 
continuity  for  an  inert  tracer  without  sources  or  sinks,  expressed  in 
flux  form  and  spherical  pressure  coordinates. 

-If  =  £  (mid  +  R  oos'i"  M  mVCOS*+R‘sh"ax"’U'  (1> 


where  m  is  the  mixing  ratio  of  the  tracer,  R  the  earth's  radius,  <b  the  geo¬ 
graphic  latitude,  X  the  longitude,  p  the  pressure,  oithe  substantial  pres¬ 
sure  derivative  following  an  individual  air  parcel,  and  u  and  v  the  wind's 
eastward  and  northward  components  respectively. 

The  local  mean  of  any  dependent  variable  in  (1)  is  defined  as  its 
average  computed  over  a  volume  characteristic  of  the  grid  distance  to  be 
used  in  the  integration.  The  value  at  any  point  within  the  region  can 
now  be  represented  by  the  sum  of  this  mean  and  a  departure  from  the  mean 

_  *  _  * 

m=m+m  u=u+u 

v  =  v  +  v*  u>  =  W  +  W*  , 
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where  the  starred  variables  are  the  departures  and  the  barred  ones  are 
the  means.  These  relations  are  then  substituted  into  (1)  and  the  entire 
relation  is  averaged  recognizing  that,  while  the  average  of  a  fluctuation 
alone  or  of  the  product  of  a  fluctuation  with  a  mean  vanish,  the  averaged 
product  of  two  fluctuations  does  not.  (See  Haltiner  &  Martin,  1957). 


=  +  t^(iBvoos4)+rdr7 


(3) 


+  (m*w*)  +  r— (m*  v*cos  <t>)  + 


ap 


R  cos  <b 


R  cos  $  ax 


(m*  u*)) 


The  terms  containing  products  of  fluctuations  are  then  approximated 
by  using  the  Prandtl  mixing-length  assumption  that  the  correlation  be¬ 
tween  velocity  and  mixing-ratio  is  proportional  to  the  gradient  of  mean 
mixing  ratio  as  given  by 


~  .  „  3m 

m*  c*j*  =  -K,  — 

i  ap 


*v*  = 


m*v 


3-  am 
R  ad 


(4) 


m*u*  =  - 


! _  am 

Rcosd  ax 


The  vertical  and  horizontal  eddy  diffusivities ,  and  are  assumed 
to  be  constant  in  both  space  and  time .  This  particular  choice  of  con¬ 
stants,  with  the  same  eddy  diffusivity  in  both  horizontal  directions  and 
without  any  off-diagonax  elements  in  the  diffusivity  tensor,  implies  that 
diffusion  is  isotropic  in  the  horizontal  plane  with  the  axis  of  least  dif¬ 
fusion  oriented  vertically  since  ^  • 

?y  substituting  (A)  for  the  eddy  correlation  terms  in  (3)  and  writing 
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the  flux  terms  back  In  advective  form  through  use  of  the  mass  continuity 
equation  for  the  local  mean  flow,  one  obtains 


3m  _  —  3m  v  £m  +  u  dm  _  „  3  m 

3t  w  3p  R  30  R  cos  0  3X  'l  .  2 

op 


K, 


R  cos  0 


(5) 


At  this  point  (5)  Is  placed  in  non-dimens tonalized  form  by  the  fol¬ 
lowing  linear  transformation  of  the  velocity  fields: 

u  =  u  u'(  v  =  v  v1,  ai  =  w  w'.  (6) 

o  o  o 

Here  the  sub-zero  values  are  constants  with  magnitudes  (u  =  v  =  20  kt, 

o  o 

0)  =  5  mb  day  J)  so  chosen  that  the  primed  variables,  whose  variations 
o 

contain  the  spatial  and  temporal  changes  of  the  velocities,  are  of  order 
one.  Since  the  equation  is  linear  in  m  no  such  scaling  needs  to  be  done 
for  that  variable  and  (5)  becomes 

v  v'  _  u  u'  .  2— 

_  am  =  u  JEL  ,  ._2 —  .  o  &D  _  K  am 

at  ow  3p  R  30  R  cos  0  3X  1  3p^ 

-  K2  (-^ - L  [cos  ] .  ,7) 

R  cos  0  30  30  R  cos"  0  3X 


For  solution  on  a  digital  computer,  the  infinitesimal  derivatives  in 
(7)  have  to  be  replaced  by  finite  differences,  thus,  converting  the  partial 
differential  equation  into  a  system  of  algebraic  equations  in  the  mixing 
ratio  and  wind  components  at  discrete  points  in  the  region  of  integration. 


This  procedure  requires  that  the  independent  variables  of  the  problem  be 
expressed  as  multiples  of  small, but  finite  increments,  sc  that  the  co¬ 
ordinates  of  a  point  (p,  0,  X,  t)  become  i,  j,kand  n;  where  i  ■  p/5p, 
j  =  {'3-dG)/6d.  k  =  X/5X  and  n  =  t/5t  and  the  values  assumed  by  (p,  <£,  X,  t) 
are  required  to  make  i,  j,  k  and  n  integers.  In  this  formalism 
m  (p,  <b,  X,t)  is  m"  and  v(p#  t)  is  vj1 

Once  such  a  coordinate  system  is  set  uft  it  becomes  possible  to 
approximate  the  time  derivative  by  a  forward  difference 


n+1  _  n 

Am.  _  mU,k  ~  mi,j  ,k 
at  at 


(8) 


the  first  order  space  derivatives  by  centered  differences,  for  example, 

(9) 


n  n  .  n 

.  —  m, , ,  .  .  -  m.  ,  ,  ,  A.  m 
3m  _  i+l,|,k  1-1, i,k _ i_ 


a  p 


26P 


25p  ' 


and  to  use  the  three  point  approximation  in  place  of  the  second  order 
space  derivatives 


2*  * 

—  m. ,  .  .  . 
d  m  _  i+1 ,  j ,  k 

2  2 
ap  (Op) 


dm"  X  A2mn 

« m .  *  ^  m ,  «  ,  <  .  ni 

l,j,k  t-l.l.k  _  _i _ 


(0p)‘ 


(10) 


The  A  operator  always  applies  at  the  point  i,  j,  k;  although  the  space  sub¬ 
scripts  are  usually  omitted  to  simplify  the  notation.  The  subscript  on  the 
operator  denotes  the  coordinate  direction  along  which  the  difference  is 
to  be  taken  and  the  superscript  the  order  of  that  difference. 

With  all  the  derivatives  replaced  by  finite  differences  and  after 
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multiplication  through  by  fit,  (7)  becomes 


m 


n+1  n  r  Wo  n  ’  L  [  vo®*  -]  A  f 

=  m  -L-sTl*  V  "LIrm  Jv  Ajm 


Upfit  uAkm  Kjfit  2  t 

L  2R6  X  J  cos  (<b  +j6ij)  +  _  (fip)^  j  ^i  m 


K-fit  A. (cos  (A  +jfii)  A.  m  ) 
+  _ i - -i  _ } _ o _ J - 

Lfi2(fi0)2  J  cos(VJ6*) 

K.fit 

r  2  _ 1 _ k 


R2(fiX)2  cos2(0+jfltf) 


(11) 


where  the  n  and  n  +  1  superscripts  denote  the  base  time  level  and  the 
newly  generated  time  level,  respectively,  and  the  l  superscript,  which 
also  denotes  time,  can  take  on  only  one  or  the  other  of  these  two  values. 

Each  of  the  factors  in  brackets  in  (14)  is  the  product  of  fit  with  a 
collection  of  constants  which,  when  combined,  has  the  units  of  recipro¬ 
cal  time.  Thu$  each  bracket  forms  the  dimensionless  ratio  of  the  time- 
step  to  a  time  characteristic  of  the  particular  process  changing  the  mixing 

ratio.  For  example,  in  the  first  term  2fip/u>  represents  the  vertical  ad- 

o 

2 

vection  time  or  in  the  fourth  term  (fip)  /K^  is  the  meridional  diffusion 
time . 

Since  each  of  these  factors  is  a  constant,  (11)  is  simplified  by  re¬ 
placing  them  with  single  quantities  to  reproduce  the  simpler  form 
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mn+1  =  mn  -  bo/ A.  m  -  b„  vi  A.n/  -  [b./cos(d  +j6<£)]i/  A.  ml 
1  1  j  3  o  k 

2  1  l 

+  b.A,m  +  [b_/cos(0  A  ,cos(<£  +j6<j)A.m 

4i  5  ojo  3 

+[bc/cos 2(<i  +jfl<a)]  A^m  f  (12) 

6  o  k 

The  pressure  increment  is  12.5  mb,  that  of  latitude  5°  and  that  of  longi¬ 
tude  10°.  With  this  information  the  first  three  advection  times  are 
evaluated  to  give: 

t  =  25  p/oj  =  5  days,  r„  =  2R5  0/v  =1.25  days, 

1  o  2  o 

r_  =  2Rfl\/u  =  2.5  days.  (13) 

3  o 

The  value  of  t 3  does  not  represent  the  actual  advection  time  enter¬ 
ing  into  the  calculation,  since  a  factor  of  cos  0  is  not  included  in  its 
evaluation.  This  factor  represents  the  decrease  in  advection  time,  as 
the  distance  corresponding  to  8X  shrinks  with  approach  to  the  pole. 

Since  the  cosine  of  the  latitude  varies  from  .819  at  35  N  to  .087  at  85 
the  actual  zonal  advection  time  changes  by  nearly  an  order  of  magni¬ 
tude  over  the  region  of  integration. 

Since  a  state  of  nearly  total  ignorance  exists  in  regard  to  the 
actual  magnitudes  of  the  eddy  diffusivities  in  the  stratosphere,  the 
diffusion  time  along  a  given  axis  is  set  to  a  simple  multiple  of  the 
advection  time  along  the  same  axis.  In  this  study,  the  multiple  (called 
!/b^)  is  chosen  to  be  eight.  This  means  that  the  parameterized  diffusion 


11-20 


can  remove  departures  from  equilibrium  1/8  ar  fast  as  advection  operating 
by  itself.  This  assumption  leads  to  the  following  values  for  the  diffusion 
times: 

=  40  days,  T5  =  5  days, 

(14) 

r  e  =  20  days . 

D 

The  difference  of  a  factor  of  four  between  r,  and  t  arises  because 

b  b 

6X  =  2  .  Since  the  horizontal  diffusion  is  assumed  to  be  isotropic, 

the  transport  must  occur  four  times  as  rapidly  over  half  the  distance. 

4  2 

The  values  of  Kj  and  K2 ,  corresponding  to  and  t  are  1  x  10  cm 
sec  1  and  4.5  x  108  cm2  sec""1  compared  to  2  x  104  and  1  x  1010  used  in 
Prabhakara's  parameterized  study.  That  author's  definition  of  diffusion 
and  the  one  used  here  are  essentially  different.  In  the  former  case  the 
diffusion  was  used  to  simulate  transfer  by  the  large-scale  eddies,  while 
in  this  study  it  has  only  to  account  for  transfer  by  fluctuations  too  small 
to  be  represented  in  synoptic-scale  data.  The  similarity  between  the 
two  vertical  diffusivities  is  unfortunate;  one  would  have  preferred  to  use 
a  vertical  diffusivity  that  differed  from  that  of  Prabhakara  by  a  factor  of 
at  least  ten  because,  in  the  parameterized  studies,  diffusion  is  called 
upon  to  provide  much  of  the  transport  that  advection  accomplishes  in 
this  model . 

The  entire  treatment  of  diffusion  is  rather  artificial,  but  since  it 
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serves  as  a  computational  smoother,  its  effect  can  not  be  minimized 
without  sacrificing  numerical  stability, 

B.  TIME  INTEGRATION 

Integration  in  time  is  accomplished  using  a  backward  corrected 
Euler  scheme  (see  Kurihara,  1965).  If  the  advective  and  diffusive  terms 
in  (12)  may  be  replaced  by  f{mn,  ojn,  vn,  u11),  this  scheme  may  be 
represented  as: 

m=m+f(m,co,v,u),  (15) 

n+1  n  ,,  l  n+1  n+1  n+1,  n 

m  =  m  +  f(m  ,  o)  ,  v  ,  u  ).  (16) 

n  A 

Given  the  field  of  m  and  the  wind  at  the  n  time  level,  m  is  computed 

l 

from  (15).  Then,  m  and  the  winds  at  the  n+1  time  level  are  used  in  (16) 
to  obtain  mn+1.  Therefore,  mn+1  is  first  predicted  using  a  forward 
difference,  and  then  corrected  with  a  backward  difference,  the  entire 
process  repeated  at  each  successive  time  level  until  m  is  known  for  the 
entire  period  of  interest. 

According  to  Kurihara  (1965),  if  this  method  is  applied  to  the  one 
dimensional  advection  equation. 


9  n  ,  ,  .  n 

x  =  x  +  bA.x 
i 

(17) 

n+1  n  ,  ,  A  l 

=  x  +  b^x  , 

(18) 

stability  will  result  for  b  i  1.  By  experiment,  it  was  discovered  that 
optimum  results  are  obtained  by  using  a  time  step  5t  =  1/48  day,  producing 
the  following  values  for  the  b's  (b_,  -  1/8): 
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b  =  4. 17  x  10~3,  b2  =  1.67  x  If)"2,  b3  =  8.33  x  10 

h  =  0.52  x  10~3,  fc  =  4.17  x  10~3,  b.  =  1 .04  x  10 
4  5  ..  b 

These  values  are  obviously  much  less  than  one,  but  the  factor  of  1/cos  d> 
makes  the  coefficient  of  the  third  term  much  larger  near  the  pole.  At 
80  N,  the  most  northerly  point  at  which  (15)  and  (16)  are  applied, 
cos  <t  =  .174  so  b„/cos  <t  =  0.047 is  still  well  within  the  region  of  stability. 

O 

Unfortunately,  a  longer  time-step  drastically  reduces  the  quality  of  the 
results.  Although  no  theoretical  investigation  is  attempted,  this 
may  be  caused  by  either  mass  imbalances  in  the  wind  data  or  by  the 
variability  of  the  space  increment  resulting  from  the  cos  <t>  effect. 

C.  WIND  DATA 

The  wind  data,  used  by  Mahlman  (1967,  1969)1  are  employed  in  the 
advective  terms  of  (11).  The  data,  which  cover  the  41  day  period  from 
15  November  until  25  December  1958,  consist  of  the  three  wind  compo¬ 
nents  and  the  fields  of  temperature  and  geopotential  height,  tabulated  at 
50  and  100  mb  for  the  region  from  40  to  80  N.  Originally, the  two  hori¬ 
zontal  wind  components,  as  well  as  the  temperature  and  height  fields, 
were  extracted  directly  from  the  Weather  Bureau  (1963)  stratospheric 
daily  chart  series  for  the  IGY,  while  the  go's  were  computed  from  the  other 
fields,  using  the  thermodynamic  equation  with  an  assumed  diabatic  heating 
rate  (Mahlman,  1967). 

During  the  early  portion  of  the  period  of  data  coverage,  the  polar 


t 

(19) 
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night  vortex  is  intensifying  with  the  circulation  dominated  by  a  cold  low 
over  northern  Asia.  Then,  on  21st  of  November  the  flow  undergoes  a  pulsa- 
tioaor  "minor  breakdown",  becoming  more  meridional  in  character  and 
remaining  disturbed  until  11  December  when  the  vortex  re-stabilizes 
with  an  important  wave  number  two  component.  During  this  last  period 
one  trough  lies  over  Asia  and  the  other  over  North  America  with  the  two 
intervening  ridges  over  the  oceans. 

With  the  exception  of  the  ten-day  period  from  5  December  until  15 
December,  the  mean  cell  in  these  data  has  rising  motion  over  the  pole 
and  descent  in  middle  latitudes  (Mahlman,  1969).  Even  during  the  time 
in  December  when  the  vertical  motion  reverses  at  the  pole,  descent  takes 
place  only  in  the  immediate  region  of  the  pole  with  rising  motion  con¬ 
tinuing  in  the  60  and  70  N  band.  This  means  that,  except  for  that  brief 
period  in  December,  the  mean  cell  operates  in  the  exact  opposite  sense 
to  that  required  by  the  Brewer-Dobson  theory. 

Unfortunately,  the  data  are  tabulated  every  24  hours  on  only  two 
pressure  surfaces,  and  for  a  somewhat  coarser  horizontal  grid.  This 
means  that  the  data  has  to  be  extended  by  extrapolation  and  interpolation 
in  both  space  and  time.  First,  using  the  thermal  wind  equation,  the  hori¬ 
zontal  winds  are  extrapolated  upward  to  12.5  mb  and  downward  to  150  mb, 
and  the  uj's  at  these  two  levels  are  assumed  to  be  zero.  This  is  done 
at  every  point  in  the  horizontal  that  carries  the  initial  data  fields  (every 
40°  of  longitude  at  80  N,  every  20°  at  70  N,  and  every  10°  of  longitude  at 
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60,  50  and  40  N).  Values  for  every  10°  longitude  at  70  N  and  80  N  are  gen¬ 
erated  by  linear  Interpolation  between  the  tabulated  values.  The  same 
process  is  applied  again  to  produce  winds  at  odd  multiples  of  five  degrees 
of  latitude.  Finally,  a  cubic  Lagrange  interpolating  polynomial  is  fitted 
to  the  data  at  each  horizontal  grid  intersection  for  the  four  levels  and 
evaluated  at  the  intermediate  pressure  levels  to  complete  the  wind  field 
for  each  day.  The  wind  for  each  time  step  is  obtained  from  the  dally 
wind  fields  by  simple  linear  interpolation  in  time. 

The  wind  data  require  still  more  treatment  before  they  are  suitable 
for  use  in  the  integration  scheme.  Due  to  inaccuracies  in  observation 
and  processing,  the  probable  errors  in  the  v  wind  component  are  at  least 
as  large  as  the  real  mean  cell.  Furthermore,  the  omegas  generated  by  the 
thermodynamic  method  are  not  necessarily  consistent  with  the  divergences 
of  the  horizontal  wind.  Mass  imbalances  In  the  large-scale  eddies,  while 
undesirable,  tend  to  cancel  each  other  when  integrated  over  the  entire 
region;  but  a  mean  cell  with  spurious  mass  sources  may  seriously  distort 
the  results  of  the  computation.  It  was  decided  to  remedy  the  situation  by 
removing  the  natural  mean  cell  and  replacing  it  with  an  artificial  time  in¬ 
dependent  one. which  exactly  satisfied  the  zonally-averaged  continuity 
equation.  The  following  transformation  accomplishes  this  end: 

u " ' .  =  u ' 

i.l,k  i,j,k 


„l  ■  l 

v  .  ,  =  V.  .  -  V.  .  +  V  , 

i,J,k  I,j,k  1,]  i , 3 


„l  - 

cu  .  =u>'  .  ,  -  a).  +  O 

i,J, k  i,j,  k  1,1  i,j 


(20) 
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Here,  the  double  primed  quantities  are  the  new  values,  the  unprimed,  the 
old,  the  tilde4  the  zonal  mear^old  values,  and  the  capitals*  the  contrived 
mean.  This  new  mean  is  usually  run  With  ascent  at  the  poles,  although 
some  experiments  were  run  with  the  cell  reversed  for  comparison.  Typi¬ 
cal  magnitudes  for  the  new  mean  were  on  the  order  of  a  knot  for  V  and 
.5  mb  day  *  for  O  .  The  new  cell  was  constant  in  time  and  no  attempt 
was  made  to  introduce  time  variations  resembling  those  of  the  natural 
mean  cell. 

While  use  of  an  artificial  mean  cell  seems  to  introduce  the  very 
element  of  arbitrariness  that  was  avoided  through  employing  observed 
winds,  it  should  be  pointed  out  that  this  approach  both  permits  experi¬ 
mentation  with  mean  cells  of  differing  character  and  allows  use  of  mean 
circulations  that  are  more  or  less  consistent  with  heat  and  momentum 
considerations . 
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III.  EXPERIMENTS 


A.  OZONE  SIMULATION 

1 .  Boundary  and  Initial  Conditions 

The  behavior  of  stratospheric  ozone  is  simulated  by  the  evolu¬ 
tion  of  an  inert  tracer  from  a  zonally  symmetric  initial  state  roughly  cor¬ 
responding  to  the  observed  zonal  mean  ozone  distribution  (Hering,  1960). 
Throughout  the  integration, the  observed  initial  values  at  12.5  and  150 
mb  are  retained  at  the  top  and  bottom  boundary  points.  This  condition 
can  be  justified, since  it  is  assumed  that  the  value  at  12.5  mb  is  con¬ 
trolled  only  by  photochemistry  and  is  time-independent  except  for  di¬ 
urnal  and  seasonal  variation  which  are  not  reproduced.  Similarly  the 
mixing-ratio  at  150  mb  is  assumed  to  represent  a  constant  equilibrium 
between  downward  transport  and  destruction  in  the  troposphere. 

Originally,  it  was  the  intention  to  parameterize  advection  into  the 
region  from  the  tropics  by  retaining  the  observed  zonal  mean  values  at 
the  southern  boundary  points ,  but  the  final  version  of  the  model  is 
permitted  to  generate  its  own  southern  boundary  condition  inasmuch  as 
the  points  at  35  N  are  all  set  to  the  generated  zonal  mean  of  those  at 
40  N. 

Near  the  pole,  the  cosine  of  the  latitude  becomes  small, causing  an 
apparent  singularity  in  the  prognostic  equation.  This  difficulty  can  be 
avoided  by  applying  a  cartesian  version  of  the  prognostic  equation  at  the 
pole;  however,  a  simpler  solution  is  suggested  by  the  treatment  of  the 
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southern  boundary,  so  the  values  at  85  N  are  set  to  the  2onal  mean  of 
those  at  80  N. 

These  boundary  conditions  imply  that  the  system  is  not  closed,  since 
boundary  fluxes  can  occur  through  all  surfaces  of  the  region  of  integration. 
A  horizontal  flux  through  the  southern  boundary  is  physically  reasonable, 
since  advection  from  tropical  regions  is  an  important  source  of  mid-  and 
high-latitude  ozone.  Similaily,  the  flux  over  the  pole  is  to  be  expected, 
although  the  boundary  scheme  used  does  not  properly  simulate  it;  a 
spurious  over-the-pole  transfer  arises  which  distorts  the  final  results. 

The  upper  boundary  condition  allows  ozone  to  be  transported  downward 
into  the  region  of  integration  but  does  not  permit  it  to  leave  through  the 
top.  The  values  at  12.5  mb  are  invariably  higher  than  those  in  the 
interior  of  the  region.  Since  the  vertical  motion  is  small  near  the  upper 
boundary,  diffusior  predominates  allowing  the  tracer  to  flow  inward  but 
never  outward.  A  similar  set  of  conditions  arises  at  the  lower  boundary; 
except  here,  the  boundary  values  are  always  lower  than  those  ir  the  in¬ 
terior  so  only  outward  fluxes  can  occur.  Thus,  although  no  source  or 
sink  terms  are  explicitly  included  in  the  model,  the  boundary  conditions 
crudely  approximate  the  behavior  of  photochemical  sources  in  the  upper 
stratosphere  with  a  sink  at  the  tropopause. 

2 .  Results 

Tigures  1  through  6  are  a  series  of  charts  portraying  the  ozone 
mixing  ratio  and  geopotential  height  on  the  50  and  100  mb  levels  at 


11-28 


Simulated  ozone  mixing  ratio  on  25  November  at  A: 


Simulated  ozone  mixing  ratio  at  100  mb  on  A:  10  December;  and  B: 


FIG.  5.  Simulated  ozone  mixing  ratio  at  100  mb  on  A:  20  December;  and  B:  25  December. 


FIG.  6.  Simulated  ozone  mixing  ratio  on  20  December  for  A:  reverse  mean  cell; 
and  B:  zero  mean  cell. 


intervals  of  5  days  during  the  period  of  integration.  Comparison  between 
the  50  mb  and  100  mb  pressure  levels  indicate  that  the  degree  of  conti¬ 
nuity  of  the  ozone  field  between  these  two  levels  is  comparable  to  that 
observed  in  the  geopotential  height  field  so  that  representation  of  more 
than  one  horizontal  section  of  the  ozone  field  is  superfluous,  and  only  the 
100  mb  surfaces  are  presented  after  2  5  November. 

These  charts  show  a  number  of  qualitative  features  that  are 
reassuringly  similar  to  the  characteristics  of  the  real  atmospheric  ozone 
field.  Dobson’s  (1924)  observations  of  ozone  maxima  in  the  troughs  are 
reproduced.  Similarly,  Bozkov's  (1968)  observations  of  sinusoidal  iso- 
pieths  of  concentration  with  ozone  maxima  in  the  standing  troughs  over 
continents  and  minima  over  the  intervening  seas  appear  in  these  fields. 

At  50  mb  and  60  N  the  normalized  correlation  coefficient  between 
to  and  m,  computed  over  the  entire  period  of  integration,  is  0,22.  This 
means  that  the  high  mixing  ratios  in  the  lows  can  be  explained  in  terms 
of  advection  from  above-  just  as  Reed  (1950)  and  Normand  (1953)  postulated. 

On  the  50  mb  surface  the  ozone  mixing  ratio  values  initially 
decrease  from  5.3  ^gm/gm  at  35  N  to  3.9  ^gm/gm  at  85  N  while  at 
100  mb  the  initial  values  increase  from  1.2  to  3 . 1  in  the  same  interval. 

As  the  integration  proceeds,  the  transport  processes  reverse  the  50  mb 
gradient  by  18  November  and  establish  a  mid-latitude  maximum  at  75  N 
by  25  November.  This  maximum  tends  to  drift  southward  appearing  most 
frequently  at  about  65  N  during  the  latter  portio'’  c:  the  experiment.  It 
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is  quite  persistent  in  this  latitude,  but  often  it  disappears  for  several 
days,  being  absent  a  total  of  19  of  the  41  simulated  days.  A  simi¬ 
lar  maximum  appears  in  the  100  mb  profile  on  2  December  but  it  is  less 
pronounced  and  not  nearly  so  persistent  in  time  (being  present  on  10  of 
41  days).  At  both  altitudes  the  maximum  zonal  mean  gradient  occurs 
south  of  the  latitude  of  maximum  mixing  ratio,  with  the  field  being 
flatter  over  the  pole.  These  results  are  consistant  with  Bozkov's  work. 
Since  they  arise  during  the  integration  and  do  not  represent  mere  per¬ 
sistence  of  the  initial  state,  they  enhance  confidence  in  the  model. 

When  time  series  of  the  zonal  mean  mixing  ratio  and  of  the 
various  processes  modifying  it  are  considered  (see  FIG.  8),  it  is  found 
that  the  horizontal  eddies  and  horizontal  diffusion  are  very  highly  cor¬ 
related  with  the  fluctuations  of  the  mixing  ratio.  The  other  processes, 
the  vertical  eddies,  the  vertical  and  horizontal  mean  cells  and  the 
vertical  diffusion  are  all  about  an  order  of  magnitude  smaller  than  these 
first  two  effects  and  did  not  produce  noticeable,  short-term  changes 
in  the  mixing  ratio. 

The  direction  of  horizontal  eddy  transport  changes  during  the 
integration.  At  the  beginning  it  is  poleward,  but  during  the  pulsation 
of  the  polar  night  vortex  it  changes  sign  and  remains  negative  for  22 
days  until  the  vortex  begins  to  re-intensify .  When  the  normalized  mixing- 
ratio  and  meridional  velocity  correlation  coefficient  is  computed  for  the 
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Simulated  zonal -mean  ozone  cross  section  in  10  gm-gm 
A:  Initial  cross  section;  B:  cross  section  for  20 
December  with  mean  ascent  at  the  pole;  C:  for  the 
same  date  with  a  zero  mean  cell;  and  D:  for  a  mean 


cell  with  descent  at  the  pole. 


entire  period,  its  value  is  +.023  at  SO  mb  and  60  N.  But  when  the 
data  are  stratified  into  poleward  and  equatorward  transport  periods,  the 
correlation  during  the  former  is  +.317  and  the  latter  -.237.  Since  the 
omega-mixing -ratio  correlation  is  +0.22,  this  indicates  that  the  dif¬ 
ference  in  the  effectiveness  between  the  vertical  and  horizontal  eddies 
is  primarily  caused  by  the  smaller  ratio  of  characteristic  velocity  to 
space  increment  for  the  former. 

Because  the  mean  cell  is  constant  in  time  and  the  zonal  mean 
gradient  is  nearly  so,  the  sign  and  magnitude  of  the  instantaneous  change 
resulting  from  the  mean  cell,  is  quite  pre  si  stent.  This  means  that,  while 
the  short-term  effect  is  small,  the  mean  cell  produces  a  marked  change 
in  the  field  when  integrated  over  41  days. 

This  effect  is  well  illustrated  by  the  differences  among  the  20 
December  100  mb  fields  (FIG.  5  and  FIG.  6).  Figure  6  shows  the  20 
December  mixing-ratio  field,  resulting  from  a  zero  mean  cell  as  well 
as  that  resulting  from  a  reversed  mean  cell.  Comparison  among  these 
three  charts  indicates  that  in  all  cases  the  features  are  very  similar  in 
shape  and  location  but  the  values  differ  as  though  the  fields  were 
translated  by  addition  of  a  constant.  This  is  very  nearly  what  happened. 
Figure  7  shows  the  zonal  mean  cross  sections  for  each  of  these  three 
cases  on  25  December  as  well  as  the  initial  zonal  mean.  Each  different 
mean  cell  changes  the  inclination  of  the  mixing  ratio  isopleths,  tending 
t.o  increase  their  slopes  when  sinking  occurs  at  the  pole  and  decrease 
them  when  rising  takes  place  at  the  polo.  This  effect  is  most  important 
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in  the  northern  portions  of  the  region  where  it  moves  the  lines  in  the 
vertical  while  toward  the  south  they  remain  at  about  the  same  altitude  in 
all  cases.  It  should  be  pointed  out  that,  while  this  effect  becomes  im¬ 
portant  only  after  long  integration  times,  it  will  eventually  attain  an 
equilibrium  with  the  other  processes  and  not  accumulate  indefinitely. 

Apparently,  the  model  seeks  to  attain  just  such  an  equilibrium 

during  this  simulation.  At  50  mb  and  60  N,  the  mixing  ratio  rises 

-7  -1 

asymptotically  to  a  value  of  about  60  10  gm/gm  (see  FIG.  8).  The  hori¬ 
zontal  eddies,  the  vertical  diffusion  and  the  horizontal  mean  all  appear 
to  contain  decaying  transients  that  largely  damp  out  by  10  December. 
Apparently  about  25  days  are  required  to  recover  frc.n  the  disequilibrium 
introduced  when  the  model  is  initiated  from  a  zonally  symetric  initial 
state . 

The  contribution  of  both  components  of  the  sub-grid-scale  dif¬ 
fusion  is  disturbingly  large.  This  phenomenon  is  dependent  upon  a 
completely  arbitrary  sub-grid-scale  parameterization  scheme  and  probably 
is  much  stronger  than  any  real  atmospheric  diffusion.  It  is  unfortunate 
that  the  diffusivities  cannot  be  reduced  without  impairing  the  model’s 
computational  properties.  The  unnaturally  strong  diffusion  may  mask 
important  advective  processes  and  spuriously  shorten  the  equilibrium 
time  as  well  as  displace  the  equilibrium  from  *ts  true  value. 
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B.  RADIOACTIVE  DEBRIS  SIMULATION 


1 .  Boundary  and  Initial  Conditions 

Clouds  of  radioactive  debris  ore  simulated  by  permitting  the 
distribution  of  tracer  to  evolve  from  an  initial  three-dimensional  Gaussian 
distribution,  centered  at  a  point  within  the  region  of  integration.  It  would 
be  more  satisfying  to  use  an  initial  state  in  which  all  the  tracer  is  con¬ 
centrated  at  a  single  grid  point,  but  in  such  expenments  an  undesirable 
computational  mode  is  excited  by  the  steep  gradients  so  that  the  results 
are  meaningless. 

The  boundary  conditions  imposed  upon  the  evolving  cloud  of 
debris  represent  a  serious  problem.  For  the  previous  experiment  with 
zonally  symmetric  ozone  fields  as  an  initial  state,  the  generated  zonal 
mean  of  points  along  the  next  interior  latitude  circle  represents  a  rea¬ 
sonable  boundary  value.  Because  the  point  source  problem  is  asymmetrical , 
constraining  the  gradient  at  the  southern  boundary  to  zero  seems  to  con¬ 
stitute  a  reasonable  boundary  condition;  but  computational  difficulties 
arise  in  this  case  when  the  wind  normal  to  the  boundary  is  strong.  The 
final  solution  to  this  problem  lies  in  a  compromise  between  the  two 
approaches.  Each  boundary  point  at  35  N  or  85  N  is  set  to  a  weighted  aver¬ 
age  of  the  zonal  mean  at  the  next  interior  parallel  and  the  value  at  the 
adjacent  point  along  that  parallel.  The  weighting  factors  are  0.7  times 
the  mean  and  0.3  times  the  adjacent  point  at  the  southern  boundary',  and 
equal  weights  are  used  at  the  northern  boundary. 
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At  the  top  and  bottom  of  the  region  there  is  little  inflow  or 


outflow  because  the  vertical  velocities  are  small.  Furthermore,  these 
regions  are  fairly  remote  from  the  centers  of  concentration.  This  means 
that  zero  gradient  at  12.5  and  at  150  mb  constitutes  a  physically  reason¬ 
able  and  numerically  stable  boundary  condition  there. 

2 .  Results 

The  evolution  of  a  cloud  injected  at  latitude  60  N  longitude  140 
W  with  initial  standard  deviation  10°  of  latitude,  20°  of  longitude  and 
2  5  mb  of  pressure  was  simulated  for  eleven  days  (See  FIGS.  9,  10).  This 
initial  position  lies  in  the  almost  purely  zonal  flow  in  the  ridge  over  the 
Northern  Pacific  so  that  the  main  mass  of  debris  is  carried  eastward 
along  the  height  contours  while  gradually  spreading  and  dispersing.  By 
23  November  the  cloud  has  migrated  about  90°  of  longitude  and  lies  in 
the  eastern  portion  of  the  Atlantic  ridge.  During  this  period  the  mixing 
ratio  at  the  last  closed  contour  decreases  from  80,  its  initial  value,  to 
20  (arbitrary  units). 

On  17  November  part  of  the  cloud  is  carried  over  the  pole  and 
then  southeastward  around  the  Siberian  low.  Between  19  and  21  November 
the  mass  rounds  the  trough  and  begins  to  move  northward  so  that  by  23 
November,  it  has  circumnavigated  the  pole  at  high  latitudes  and  re-merged 
with  the  main  cloud. 

In  the  zonal  mean,  the  center  of  the  initial  distribution  lay  at 
60  N  and  a  pressure  altitude  of  75  mb  (See  FIG.  11).  As  the  cloud  evolves. 
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FIG.  10.  Simulated  radioactive  debris  at  100  mb  on  A:  22  November;  and  B:  24  November. 


FI G.  11. 


Zonal -mean  simulated  radioactive  debris  in  arbitrary 
units  on  A:  18  November;  B:  20  November;  C-  22  7 
November;  0:  24  November;  and  D:  26  November 
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the  center  of  maximum  mixing  ratio  gains  altitude  and  move:;  toward  the 
south,  at  the  same  time  spreading  laterally  in  such  a  manner  that  its 
axis  slopes  downward  toward  the  pole.  By  23  November,  the  zonal  mean 
bears  a  striking  resemblance  to  both  the  results  of  parameterized  studies 
and  to  observations  (Gudiksen  et  al.  1969),  showing  that  this  model 
can,  through  advection  and  diffusion  alone,  quickly  shape  any  reasonable 
distribution  of  inert  tracer  into  a  fair  reproduction  of  observed  mixing 
ratio  fields.  In  fact,  this  process  took  place  much  more  quickly  than  it 
would  have  in  the  real  atmosphere  because  of  the  large  standard  deviation 
in  the  initial  distribution. 

No  detailed  analysis  is  made  of  the  principal  agencies  responsible 
for  the  tracer's  evolution,  but  the  horizontal  eddies  and  diffusion  com¬ 
bined  are  about  an  order  of  magnitude  stronger  than  the  sum  of  the  other 
four  agencies.  Therefore,  as  in  the  ozone  experiments,  the  sign  of  the 
tendencies,  produced  oy  the  mean  cell  terms,  remain  consistent  in  time  so 
that  over  periods  longer  than  eleven  days  their  effect  becomes  significant. 
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IV.  CONCLUSIONS 


This  model  successfully  simulates  many  observed  qualitative  features 
of  stratospheric  transports  of  both  ozone  and  radioactive  debris.  Apart 
from  some  numerical  difficulties  which  result  from  mass  imbalance  in  the 
wind  field  and  from  inadequate  treatment  of  the  singularity  at  the  pole, 
no  obstacles  exist  to  prevent  the  use  of  this  model  to  predict  the  evolu¬ 
tion  of  stratospheric  tracers. 

If  the  erroneous  exaggeration  of  the  parameterized  diffusion  may  be 
neglected,  the  instantaneous  derivative  is  primarily  determined  by  the 
horizontal  eddies.  Although  the  effect  of  the  vertical  eddies  is  about 
an  order  of  magnitude  weaker  than  that  of  the  horizontal  eddies,  it  cannot 
be  neglected  in  the  zonal  mean.  The  two  components  of  the  mean  cell 
are  still  weaker  than  the  vertical  eddies.  However,  over  long  periods  of 
time,  they  may  assume  greater  importance  because  they  do  not  change 
sign  or  fluctuate  greatly  in  magnitude. 
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Abstract 


The  vertical  distribution  of  ozone  was  observed  by  the  Air  Force 
Cambridge  Research  Laboratories  North  American  ozonesonde  network  from 

29  April  1963  to  10  May  1963  on  a  daily  basis.  This  observed  data, 
supplemented  by  twenty-four  hour  trajectory  calculations,  was  used  to 
prepare  distributions  of  ozone  mixing  ratio  on  the  200,  100,  50  and 

30  mb  surfaces  .  The  ozone  distributions  were  evaluated  for  the  calcula¬ 
tion  of  ozone  transport  by  mean  and  eddy  motions  in  the  lower  stratosphere 
The  results  of  these  calculations  show  the  magnitude  of  the  horizontal 
and  vertical  transport  of  ozone  into  and  out  of  the  30/50,  50/100,  and 
100/200  mb  volume  over  North  America.  The  mean  vertical  transport, 
supported  by  the  vertical  eddy  transport,  is  combined  to  move  the  ozone 
primarily  downward  through  the  stratosphere.  At  the  tropopause,  the 
vertical  eddies  transport  the  ozone  to  the  troposphere  for  eventual 
destruction. 
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I.  INTRODUCTION 

Ozone  is  an  important  but  not  well-documented  atmospheric  gas.  Its 
importance  to  meteorologists  is  twofold.  First,  ozone  can  be  used  as  a 
quasi-conservative  atmospheric  trace  substance  below  30  kilometers; 
secondly,  as  an  absorber  of  ultraviolet  light,  ozone  affects  the  heat  balance 
through  the  diabatic  heating  term  of  the  first  law  of  thermodynamics.  Although 
the  exchange  of  energy  by  ozone  is  worthy  of  investigation,  ozone,  as  an 
atmospheric  tracer,  is  the  subject  of  this  study.  Before  proceeding,  a 
review  of  the  significant  works  of  earlier  investigators  is  desirable. 

The  first  major  investigation  of  the  measurement  and  distribution  of 
atmospheric  ozone  was  carried  out  under  the  direction  of  G.  M.  B.  Dobson 
during  a  five-year  period  from  1924  to  1929.  Dobson  and  his  associates 
reported  the  findings  of  their  study  in  a  series  of  four  articles  in  the 
Proceedings  of  the  Royal  Society  of  London  (1926,  1927,  1929,  1930). 

Based  on  the  work  of  previous  investigators,  Dobson  developed  a  spectroscopic 
instrument  for  use  on  the  earth's  surface  that  measures  the  total  amount  of 
atmospheric  ozone  in  a  vertical  column.  The  results  of  daily  measurements 
of  total  ozone  at  Oxford  during  1924  and  1925  were  compared  with  surface 
pressure.  It  was  found  that  surface  pressure  had  a  high  negative  correlation, 
on  the  order  of  -.5,  with  total  ozone  measurements .  On  an  annual  basis,  a 
maximum  of  total  ozone  was  observed  in  the  spring  with  a  corresponding 
minimum  in  the  fall.  Comparisons  of  total  ozone  were  made  with  atmospheric 
pressure  at  nine,  twelve,  and  fourteen  kilometers.  A  high  negative  correlation, 
on  the  order  of  -.o,  was  observed. 


Ill— 9 


In  1926  instruments  were  sent  to  various  locations  in  northwestern 
Europe.  The  observations  produced  similar  conclusions  as  the  Oxford  findings 
of  1925.  Using  the  results  of  1926-1927,  a  complete  description  of  the  ozone 
conditions  relative  to  cyclones  and  anticyclones,  was  presented.  In  general, 
high  total  ozone  was  measured  in  the  area  of  northerly  wind  flow  and  low  total 
ozone  in  the  sectors  associated  with  southerly  wind  flow.  Highest  total  ozone 
values  were  found  to  the  southwest  of  the  surface  low  pressure  centers. 

During  1928  and  1929  the  spectrographs  were  redistributed,  so  that 
observations  were  taken  at  Arosa  (Switzerland),  Table  Mountain  (California), 
Helwan  (Egypt),  Kodaikanal  (India),  Christchurch  (New  Zealand),  Oxford 
(England),  and  Montezuma  (Chile).  The  analysis  of  the  measurements  showed 
that  the  midlatitude  stations  had  similar  observations  as  discussed  for  north¬ 
western  Europe.  The  tropical  stations  observed  no  significant  pressure  or 
total  ozone  variations. 

Several  years  later,  using  the  observations  obtained  during  the  period 
1928-1932  in  northwestern  Europe,  Meetham  (1937)  found  a  high  positive 
correlation  between  total  ozone  and  temperature  at  a  height  of  12,  15,  and 
18  kilometers.  The  magnitude  of  the  correlation  coefficients  was  on  the 
order  of  +.5.  He  also  found  a  high  negative  correlation  between  total  ozone 
and  pressure  at  9,  12,  15,  and  18  kilometers  altitude  and  between  total  ozone 
and  the  height  of  the  tropopause.  Although  values  vary  somewhat,  the 
pressure  correlation  coefficients  have  a  magnitude  of  -.4,  while  the  corre¬ 
sponding  tropopause  figures  are  -  .5. 

Godson  (1960)  has  summarized  Northern  Hemisphere  total  ozone  obser¬ 
vations  OV'T  many  years  by  month  and  latitude.  His  results  show  a  strong 
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middle  and  north  latitude  minimum  during  the  month  of  October.  Relatively 
constant  values  over  the  seasons  are  shown  in  the  tropics  with  increasing 
values  northward.  The  maximum  value  is  reported  above  70N  during  the 

■j 

month  of  March.  The  steepest  latitudinal  gradient  occurs  across  middle  lati¬ 
tudes  during  late  winter  and  early  spring.  The  flattest  gradient  occurs  during 
late  summer  and  early  autumn.  This  elaborates  upon  and  confirms  earlier 
and  limited  observations.  Godson's  summary  is  generally  accepted  as  a 
reasonable  description  of  the  mean  distribution  of  total  ozone. 

The  vertical  distribution  of  ozone  has  been  difficult  to  obtain.  This  is 
attributed  to  the  lack  of  a  reliably  accurate/  inexpensive,  and  expendable 
instrument.  A  summation  of  the  various  methods  of  measuring  the  vertical 
distribution  of  ozone  is  presented  by  Craig  (1965),  and  therefore  is  not 
presented  here .  The  most  widely  used  technique  has  been  the  Umkehr  method 
devised  by  F.  P.  W.  Gotz  in  1931  (Craig,  1965).  Using  photometery,  this 
method  breaks  the  atmosphere  into  five  or  nine  layers.  For  each  layer  a 
single  ozone  value  is  computed.  Of  more  recent  use  are  the  Regener  chemilum¬ 
inescent  ozonesonde  and  the  Brewer-Mast  ozonesonde.  Both  of  thes  instru¬ 
ments  are  capable  of  being  used  on  radiosondes  and  are  able  to  give  a  much 
finer  vertical  ozone  structure  than  the  Umkehr  method . 

Investigators  (see,  for  example,  Gotz,  1951;  Ramanathan  and  Kulkarni, 
1960;  Kulkarni,  1962;  Mateer  and  Godson,  1960;  Craig,  1960,  1965;  Hering 
and  Borden,  1965;  or  Bojkov,  1965,  1969)  agree  on  the  general  shape  of  the 
vertical  ozone  profile.  Low,  relatively  constant,  mixing  ratios  of  ozone  are 
found  in  the  troposphere  with  a  very  sharp  increase  at  the  tropopause.  The 
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maximum  value  of  ozone  mixing  ratio  is  most  often  found  in  the  15  to  30 
kilometer  layer  with  the  altitude  of  the  maximum  decreasing  poleward. 

Above  this  maximum  the  ozone  mixing  ratio  slowly  decreases  or  remains 
relatively  constant  with  altitude.  Maximum  hemispheric  values  are  observed 
near  the  30  kilometer  altitude  of  the  equatorial  latitudes.  In  the  stratosphere 
isolines  of  ozone  distribution  normally  slope  poleward  and  downward  from 
the  tropical  regions .  Large  daily  fluctuations  in  ozone  content  are  observed 
in  the  region  between  the  tropopause  and  the  primary  ozone  maximum.  The 
fluctuations  are  weakest  in  tropical  regions  and  strongest  in  middle  and 
north  latitudes. 

It  is  worthwhile  at  this  point  to  consider  the  expected  distribution  of 
ozone,  assuming  photochemical  equilibrium .  Gotz  (1951),  Craig  (1950,  1965) 
Prabhakara  (1963),  and  Brewer  and  Wilson  (1968)  have  discussed  the  origin 
and  details  of  photochemical  theory  with  respect  to  atmospheric  ozone.  There 
fore,  photochemical  theory  will  not  be  reviewed  here.  From  the  work  of  the 
previously  noted  authors,  two  salient  points  are  noted: 

1 .  Photochemical  equilibrium  theory  correctly  predicts  the  observed 
amount  of  ozone  in  the  atmosphere  above  approximately  30  kilometers. 
Predicted  by  this  theory  is  the  equatorial  maximum  region  and  decreasing 
values  poleward  and  with  altitude.  Below  this  region,  photochemical  theory 
does  not  agree  with  observed  distributions. 

2.  Production  and  recombination  can  be  considered  small  below  the  30 
kilometer  maximum  observed  values.  Therefore,  in  this  region  ozone  must  be 
quasi-con  orvative  and  can  be  used  as  an  atmospheric  tracer. 
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Since  the  observed  distribution  and  photochemical  theory  disagree,  the 
discrepancy  must  be  due  to  transport  processes.  Transport  mechanisms  must 
move  the  ozone  from  high  equatorial  source  regions  to  the  lower  stratosphere 

j 

in  middle  and  high  latitudes.  From  the  lower  stratosphere  ozone  is  passed 
to  the  troposphere  for  eventual  destruction  with  atmospheric  constituents  or 
the  earth's  surface  (Craig,  1950,  1965).  Craig  further  mentions  that  the 
process  would  have  to  be  "most  effective  in  middle  and  high  latitudes  in 
winter,  when  the  total  amount  of  ozone  increases  most  rapidly"  (Craig, 

1950). 

Prabhakara  (1963)  performed  a  study  of  a  mathematical  model  of  ozone 
distribution  based  on  photochemical  equilibrium  below  41  kilometers,  and  the 
effect  of  transport  processes  upon  this  distribution.  His  conclusion,  that 
"photochemical  theory  by  itself  cannot  explain  the  latitudinal  and  seasonal 
variation  of  ozone,  '  lends  further  support  to  ideas  advanced  by  Craig  (1950, 
1965). 

In  the  past,  sufficient  ozone  measurements  on  a  hemispheric  and  vertical 
scale  were  unavailable.  This  prevented  direct  calculations  of  the  transport 
of  ozone  by  dynamic  processes  in  the  stratosphere.  Newell  (1961,  1962,  1964) 
advanced  the  idea,  that  since  the  total  ozone  measured  at  the  surface  was 
highly  correlated  with  the  ozone  amount  in  the  12  to  24  kilometer  layer,  then 
a  measurement  of  the  flux  of  total  ozone  would  be  indicative  of  the  flux  in  the 
lower  stratosphere.  Using  100  and  50  mb  winds  and  total  ozone  measurements 
during  the  International  Geophysical  Year,  Newell  computed  the  mean  and  eddy 
fluxes  of  ozone  in  this  layer.  Newell  concluded  that  the  mean  meridional 
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circulation  and  standing  eddies  did  not  serve  as  effective  transport  mech¬ 
anisms.  Thus,  he  argued  that  most  ozone  transport  was  attributed  to  quasi¬ 
horizontal  transient  eddies. 

To  explain  the  observed  change  in  total  ozone  measured  at  the  surface, 
a  mechanism  was  proposed  by  Reed  (1950)  and  Normand  (1953).  Although 
proposed  separately,  this  theory  is  usually  called  the  "Reed-Normand  effect." 

On  a  synoptic  scale  the  mechanism  combines  the  horizontal  advection  and 
vertical  motion  to  explain  the  variation  of  total  ozone  with  the  synoptic 
features  of  the  upper  troposphere  and  lower  stratosphere.  Basically  the 
argument  is  as  follows. 

Air  blows  through  the  ridge  and  trough  patterns  of  the  height  field.  The 
temperature  pattern,  associated  with  the  height  field  in  the  lower  strato¬ 
sphere,  is  fairly  constant  with  warm  areas  in  the  base  of  troughs  and  cold 
areas  in  the  ridges.  Therefore,  the  quasi-stationary  temperature  patterns 
must  be  maintained  by  adiabatic  warming  and  cooling  of  air  parcels  moving 
through  them.  The  subsidence  associated  with  the  warming,  brings  down  ozone- 
enriched  air  from  above.  The  lifting,  associated  with  the  cooling,  transports 
ozone-poor  air  from  lower  altitudes .  Since  this  region  is  below  the  level  of 
the  observed  equatorial  maximum,  higher  ozone  values  are  brought  southward 
to  the  trough  and  lower  values  are  brought  northward  from  equatorial  regions 
to  the  ridge . 

The  basic  qualitative  theories  of  movement  and  approximate  distribution 
of  ozone  in  the  lower  stratosphere  have  been  proposed.  Quantitative  conformation 
of  these  theories  has  been  lacking,  due  to  a  minimum  of  information  with 
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respect  to  the  horizontal  and  vertical  distribution  of  ozone  in  the  lower 
stratosphere  (and  elsewhere).  With  the  challenge  thus  defined,  this  study 
was  undertaken  to  compare  earlier  theories  with  a  quantitative  description  of 
ozone  movements . 

As  a  first  step,  the  distribution  of  ozone  mixing  ratio  on  standard  pressure 
surfaces  in  the  lower  stratosphere  was  constructed  from  the  observations 
taken  by  the  Air  Force  Cambridge  Research  Laboratories  ozonesonde  network 
during  29  April  1963  to  10  May  1963.  It  follows  that  ozone  motions  can  be 
evaluated,  using  basic  meteorological  data. 

To  accomplish  this  evaluation,  -  computational  scheme  was  proposed  by 
Professor  J.  D.  Mahlman  of  the  Deptutment  of  Meteorology,  Naval  Postgraduate 
School.  This  scheme  puts  the  c  acuity  equation  for  ozone  (Haltiner  and 
Martin,  1957)  into  a  form  that  enables  computation  of  the  mean  and  eddy  flux 
of  ozone  mixing  ratio  into  a  volume  bounded  by  two  pressure  surfaces.  The 
intent  of  the  study  is  to  perform  quantitative  calculations  of  ozone  transport, 
due  to  various  physical  processes  based  upon  this  formulation  and  assess 
their  significance  in  terms  of  previous  estimates. 
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11  •  -PREPARATION  OF  D/.TA  AND  THE  COMPUTATIONAL  SCHEME 
A.  BACKGROUND 

In  1963  an  attempt  was  made  to  gather  a  more  complete  representation 
of  the  distribution  of  ozone  on  a  synoptic  scale.  The  Air  Force  Cambridge 
Research  Laboratories  (AFCRL)  operated  an  eleven-station  ozonesonde  network 
in  North  America.  Daily  observations,  obtained  by  this  network  between 
29  April  1963  and  10  May  1963,  provide  the  basic  data  for  this  study  (Hering 
and  Borden,  1964).  This  period  is  just  after  the  observed  spring  maximum  of 
total  ozone.  Wind  and  temperature  fields  at  standard  levels  were  available 
from  the  National  Weather  Records  Center,  Asheville,  North  Carolina  and 
from  the  Free  University  of  Berlin  {1963,  1966). 

B.  PHYSICAL  LIMITATIONS 

Considering  the  geographical  locations  of  the  ozonesonde  stations  on  the 
North  American  continent,  an  upper  boundary  of  24  kilometers  defined  by 
sufficient  ozonesonde  observations,  and  a  lower  boundary  defined  by  the 
tropopause,  the  volume  in  which  calculations  can  be  made  is  limited. 

Four  pressure  levels  -  200,  100,  50  and  30  mb  -  were  chosen  on  the  basis 
of  readily  available  wind  and  temperature  data.  The  selection  of  the  grid 
boundaries  at  30N-140W,  75N-140W,  75N-60W,  and  30N-60W  minimizes  the 
Atlantic  Ocean,  Pacific  Ocean,  and  Northern  Canada  data-sparse  areas. 

A  latitude-longitude  oriented  grid  was  used  in  order  to  measure  fluxes 
across  boundaries  of  meridians  and  parallels.  The  grid  interval  on  a  pressure 
surface  is  fve  degrees,  which  is  sufficient  to  represent  synoptic  scale  motions 
in  the  stratosphere. 
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For  the  four  isobaric  surfaces  involved  the  distance  between  them  is 
on  the  order  of  four  kilometers.  Considering  the  vertical  distribution  of  ozone, 
this  interval  is  larger  than  desired  and  forces  rough  approximations  in  com¬ 
puting  vertical  derivatives. 

The  time  interval  was  chosen  as  one  day  corresponding  to  the  frequency 
of  ozonesonde  observations.  Considering  the  expected  movement  of  ozone 
in  this  region  and  the  synoptic  scale  of  the  analysis,  this  interval  is 
reasonable . 

C.  BASIC  ATMOSPHERIC  PARAMETERS 

Hand  analyses  of  streamline,  isotach,  and  isotherm  patterns  were  prepared 
at  each  pressure  surface  for  the  time  period  of  the  calculations.  The  daily 
2C0  and  100  mb  1200  GMT  Weather  Bureau  synoptic  analysis  and  the  daily 
50  and  30  mb  0000  GMT  synoptic  analysis  from  the  Free  University  of  Berlin 
were  used  as  a  plotted  data  source. 

As  shown  in  Appendix  A,  the  streamline  analysis  describe  typical  patterns 


^  r  ‘  'PU  ^  I  A»  /*»  •»  I  A*  1  «  j  II  «  —  k-»-  V-  — —  -  .  -  ~  k  -  —  —  /lL  —  V 

*-■*■*..  •  ihv  IV'*  'o*  *  W*wl  »  U  LliUC 


weaken  with  altitude.  At  the  30  mb  surface  the  weakening  winter  circumpolar 
vortex  is  observed  with  weak  anticyclonic  centers  forming  and  dissipating  at 
the  southern  latitudes.  The  ten-day  period  of  this  investigation's  major 
calculation  is  during  the  spring  reversal  of  the  upper  stratospheric  circulation 
from  winter  westerlies  to  summer  easterlies. 

Since  the  ozonesonde  observations  were  taken  at  1200  GMT,  an  approx¬ 


imation  to  the  parameters  of  the  atmosphere  at  the  50  and  30  mb  surfaces  was 
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made  as  follows.  The  zonal  and  meridional  components  of  the  wind  field  and 
the  temperature  field  were  linearly  interpolated  for  1200  GMT  values  from  the 
0000  GMT  fields. 

Vertical  motions  were  computed  by  solving  the  first  law  of  thermodynamics 
in  the  form 


u; 


aT 

at 


+  \V0»  VT  -  ~ 
z  c 

_ E 


dH 

dt 


a_  3T 
cp  *  3p 


(1) 


The  diabatic  heating  term  was  evaluated,  using  the  net  heating  rates,  as 
proposed  by  Kennedy  (1964)  for  the  spring  period.  Unrepresentative  changes 
in  vertical  motions  between  grid  points  and  excessive  magnitudes  at  certain 
grid  points  are  the  result  of  computational  inaccuracies  attributed  to  the 
assumption  of  a  linear  change  of  wind  and  temperature  across  a  horizontal 
grid  interval.  To  avoid  these  unrepresentative  and  dynamically  untenable 
vertical  motions,  a  smoothing  scheme  was  used  as  follows: 


"i.j  =  +  +  +  asVi.j  + 


(2) 


To  aid  in  establishing  the  state  of  the  atmosphere  for  the  area  of  interest 
during  the  period  of  the  calculations,  ten  day  latitudinal  averages  by  pressure 
surface  were  prepared  for  the  temperature  field  and  the  zonal,  meridional, 
and  vertical  components  of  the  wind  fields.  These  values  are  presented  in 
meridional  cross  section  format  as  figures  1-4. 

The  values  in  figures  1-4  are  in  general  agreement  with  results  presented 
by  other  authors  (e.g,  Craig,  1965).  The  temperature  field  (fig.  1)  is  typical 
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FIC.  2.  Ten-dav  Zonal  Hem  u-component  of  Wind  (Knots)  for  the 
30  April  to  9  May  1963.  Computed  between  longitude  60W  and  140W 
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for  the  lower  stratosphere  at  this  time  of  year.  Cold  temperatures  at  southern 
latitudes  give  way  to  the  warm  temperatures  associated  with  the  middle  and 
northern  latitude  troughs.  This  latitudinal  temperature  gradient  weakens 
with  altitude.  The  u-component  of  the  wind  (fig.  2)  supports  the  earlier 
observation  that  the  winter  circumpolar  vortex  is  weakening  and  being 
replaced  from  above  by  summer  easterlies  as  can  be  seen  at  30  mb  and  35N. 
The  mid-latitude  wind  maximum,  associated  with  the  mean  polar  front  jet 
position,  is  seen  at  200  mb  and  45N.  The  v-component  of  the  wind  (fig.  3) 
shows  little  motion  at  the  30  and  50  mb  levels.  At  100  mb  strong,  southerly 
movement  north  of  65N  and  strong,  northerly  movement  south  of  SON  is 
observed.  At  200  mb  strong,  southerly  motions  are  seen  north  of  65N.  The 
vertical  motions  (fig.  4}  show  a  general  subsidence  over  the  volume  of 
calculations  except  for  a  weak  rising  area  at  50  mb  and  strong  upward  motions 
between  45N  and  50N  at  200  mb. 

Most  investigators  presenting  similar  results  use  time  scales  of  months 
or  seasons.  This  study  is  limited  to  ten  days,  and  hence  the  values  are  not 
directly  comparable.  In  addition,  to  properly  discuss  the  various  details  of 
these  mean  atmospheric  parameters,  a  separate  and  more  thorough  treatment 
is  required  than  can  be  presented  here.  Therefore,  further  discussion  of 
figures  1  -  4  is  not  presented . 

D.  OZONE  DATA  AND  ANALYSIS 

1 .  The  AFCRL  Ozonesonde  Network 

For  most  of  its  duration  the  AFCRL  ozonesonde  network  made  weekly 
observations  on  Wednesday  except  for  selected  periods  when  more  frequent 
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observations  were  made.  During  one  such  period  from  29  April  1963  to 
10  May  1963,  daily  observations  were  taken  at  1200  GMT.  Participating 
stations  are  listed  in  Table  1. 


OZONESONDE  NETWORK 
(after  Hering  and  Borden,  1964) 


STATION 

I3SK9I 

■E3K5H 

Albrook  Field,  Canal  Zone  (AWS) 

9.0 

79.6 

Colorado  State  University,  Fort  Collins 

40.6 

105.1 

Eielson  AFB,  Fairbanks,  Alaska  (AWS) 

64.8 

147.9 

Florida  State  University,  Tallahassee 

30.4 

84.3 

Fort  Churchill,  Manitoba  (Canadian  Met.  Br.) 

58.8 

94.1 

Goose  Bay,  Labrador  (Canadian  Met.  Br.) 

53.3 

60.4 

L.  G.  Hanscom  Field,  Bedford,  Massachusetts 

42.5 

71.3 

Thule  AFB,  Greenland  (AWS) 

76.5 

68.8 

University  of  Mew  Mexico,  Albuquerque 

35.0 

106,6 

University  of  Washington,  Seattle 

47.4 

122.3 

University  of  Wisconsin,  Madison  ^ 

43.1 

■  ■  ■ 

89.4 

TABLE  I 


The  instrument,  used  in  determining  the  amount  of  ozone  in  the  atmo¬ 
sphere,  was  the  dry  chemiluminescent  ozonesonde,  developed  by  Regener 
(1960).  The  ozonesonde  was  attached  to  a  sounding  balloon  with  a  standard 
radiosonde  package.  The  instrument  measures  the  amount  of  o zone  by 
drawing  ozone-laden  air  across  a  disk  coated  with  a  chemical  that  produces 


III— 2  4 


a  luminescent  reaction  with  ozone.  The  intensity  of  the  reaction  is  measured 
by  a  photomultiplier  tube  and  in  turn  is  telemetered  back  to  ground  using 
standard  radiosonde  equipment  (Regener,  I960,  1964).  Samples  were  taken 
every  15  seconds  for  a  vertical  resolution  of  approximately  250  feet  (Hering 
and  Borden,  1964). 

2 .  Reliability  of  Measured  Chemiluminescent  Ozonesonde  Data 
In  view  of  the  uncertainties  inherent  wiiu  the  measurement  and 
calibration  of  the  Regener  ozonesonde,  it  is  worthwhile  to  note  the  major 
corrections  applied  to  the  soundings  and  to  note  the  results  of  several  inter- 
comparison  tests  with  other  ozonesondes. 

The  disk,  coated  with  the  chemiluminescent  substance,  is  sensitive 
to  light,  moisture,  and  ozone.  Therefore,  as  a  check,  calibration  shortly 
before  launch  Is  required  (Regener,  1964).  Hering  and  Borden  (1964)  report 
that  achievement  of  reliable  preflight  calibration  was  difficult.  This  in  turn 
required  a  correction  to  the  telemetered  data .  The  correction  was  achieved 


by  adjusting  the  integrated  observed  ozone  density  to  a  simultaneously  observed 
total  ozone  measurement,  using  a  spectrophotometer.  The  ozonagrams  pub¬ 
lished  by  Hering  and  Borden  (1964)  have  this  correction  applied. 

Hering  and  Dutsch  (1965)  compared  the  results  of  simultaneous 


soundings  of  the  Rogener  cheruilutni.iiesv.eui  ozonesonde  and  uiie  bubbler-type 


electrochemical  ozonesonde.  Their  conclusion  is  that  "ozonesondes  of  the 


types  tested  are  capable  of  providing  reliable  high-resolution  measurements 


of  the  vertical  ozone  distribution 


H 


More  recently  Komhyr,  Grass,  and  Proulx  (1968)  made  intercomparison 
tests  in  various  combinations  between  the  Regener  chemiluminescent 
ozonesondes,  carbon-iodine  ozonesondes,  and  Brewer- Mast  ozonesondes. 

The  result  of  the  comparisons  after  corrections  was,  that  the  chemiluminescent 
ozonesonda  compared  with  the  other  ozonesondes  in  the  stratosphere  and  gave 
"relatively  reliable"  results. 

Hering  and  Borden  (1964)  considered  the  data,  obtained  by  the  ozone- 
sonde  network,  as  "provisional."  Considering  the  results  of  the  two  inter¬ 
comparison  tests  and  the  observation  that  other  investigators  (Berggren  and 
Labitzke ,  1966,  1968;  Breiland,  1964,  1965,  1967,  1968;  Craig,  Deluisi, 
and  Sticksel,  1967;  Hering,  1966;  Hering,  Touart,  and  Borden,  1967;  and 
Sticksel ,  1966)  have  published  studies  using  ihe  corrected  ozonesonde  data, 
the  reliability  of  the  ozonesonde  soundings  as  published  by  Hering  and 
Borden  (1964)  is  accepted. 

3 .  Preparation  of  Ozone  Data  for  Analysis 

Although  the  ozone  soundings  were  accepted,  there  are  several 
features  of  the  soundings  that  had  to  be  modified  to  be  consistent  with  the 
scale  of  the  intended  calculations.  On  the  ozonagram  thin  layers  with  large 
unrepresentative  gradients  and  layers  in  which  the  ozone  values  fluctuated 
excessively  and  with  apparent  randomness  were  considered  not  representative 
of  the  synoptic  scale  distribution.  These  features  of  the  vertical  distribution 
may  either  be  a  failure  of  the  instrument  or  a  smaller-scale  ozone  distribution 
which  is  presently  below  our  ability  to  identify.  To  overcome  these  apparent 
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inconsistencies,  a  subjective  smoothing  of  the  individual  soundings  was 
completed  as  the  first  step  in  analysis.  This  smoothing  is  not  to  be  over¬ 
emphasized  since  extreme  care  was  taken  not  to  destroy  the  characteristics 
of  the  soundings.  The  severest  smoothing  was  on  the  order  of  .3  kilometers 
on  the  height  scale  with  95  percent  of  the  ozonagram  values  being  unaffected 
by  the  smoothing  process. 

With  the  smoothed  values  obtained,  a  time  series  for  the  twelve-day 
period  wa~  prepared  for  each  station.  Linearly  interpolated  values  were 
calculated  for  stations  missing  an  observation,  but  having  an  observation  24 
hours  prior  to  and  after  the  missing  observation.  Linear  interpolations  were 
also  made  in  the  case  of  the  Fairbanks  station  where  observations  were 
reported  at  0000  GMT. 

To  supplement  the  small  number  of  observations  for  the  area  under 
consideration,  hand  trajectories  were  calculated  for  the  24-hour  periods  prior 
to  and  following  an  observation.  This  technique  used  12-hour  streamline  and 
isotach  patterns  to  obtain  the  position  of  the  parcel  24  hours  before  and  after 
the  observation  time.  Under  conservation  of  ozone  mixing  ratio  ),  these 
new  points  then  possess  a  value  of  x  corresponding  to  the  value  measured  at 
the  original  point.  An  estimate  of  the  vertical  motion  during  the  trajectory  was 
made,  using  an  approximation  derived  from  the  Eulerian  expansion  of  the  total 
derivative  of  pressure  with  respect  to  time,  the  hydrostatic  equation,  and  the 
first  law  of  thermodynamics.  With  this  estimate  of  the  vertical  motion  applied 
to  the  initial  ozonagram  sounding,  a  corrected  trajectory  value  was  obtained. 
The  trajectory  calculations  are  justifiable,  assuming  the  conservative  property 
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of  ozone  mixing  ratio.  As  linearly  interpolated  values  were  not  used  for 
trajectory  computations,  this  calculation  increased  the  number  of  data  points 
between  two  and  three  times  the  original  number  of  ozone  data  points 
available  for  analysis  on  a  pressute  surface. 

4.  Hand  Analysis  Techniques 

The  manner  in  which  the  data  was  to  be  analyzed  presented  somewhat 
of  a  problem.  Synoptic  models  of  ozone  distribution  on  a  pressure  surface 
have  not  been  published  to  the  author's  knowledge — with  the  exception  of  the 
six  analyses  presented  by  Berggren  and  Labitzke  (1968).  Other  investigators 
have  published  information  about  the  synoptic  distribution  of  ozone  but  none 
have  presented  distributions  on  a  pressure  surface.  Pittock  (1969)  has 
investigated  the  synoptic  climatology  of  ozone,  using  statistical  methods. 

His  "preliminary  results, "  using  546  Brewer-Mast  ozonesonde  soundings,  are 
primarily  in  tabular  format,  classed  into  general  synoptic  conditions  by 
season.  Breiland  (1964,  1965,  1967,  1968)  studied  the  distribution  of  ozone 
on  a  sub-synoptic  scale,  primarily  using  vertical  cross  sections,  time 
sections,  and  individual  AFCRL  network  soundings.  Sticksel  (1966)  analyzed 
the  vertical  distribution  of  ozone  with  data  collected  at  the  Tallahassee  station 
of  the  AFCRL  network.  Briggs  and  Roach  (1963)  and  Danielsen  (1968)  have 
investigated  the  ozone  distribution,  along  with  other  parameters,  near  the  jet 
stream  and  the  tropopause,  using  aircraft  observations. 

The  first  attempt  to  analyze  the  plotted  data  points  was  inconclusive. 
Merely  constructing  isolines  of  ozone  mixing  ratio  from  plotted  data  points 
left  too  much  to  the  imagination.  This  resulted  in  inconsistent  analyses  when 
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compared  to  the  work  of  the  previously  mentioned  authors .  These  analyses 
were  also  in  consistent  in  time. 

To  overcome  this  problem  the  appropriate  streamline  analysis  was  put 
under  the  plotted  ozone  data.  With  the  streamlines  to  guide  the  analyst  in 
sparse  data  areas,  the  desired  consistency  was  achieved.  To  elaborate 
further,  the  criteria  for  placement  of  isolines  of  ozone  mixing  ratio  are 
listed  below  in  order  of  priority: 

a)  observed  data  at  map  time; 

b)  data  from  trajectory  calculations  as  corrected  for  vertical 
motion; 

c)  streamline  pattern; 

d)  vertical  motion  patterns. 

Justification  for  this  priority  is  obvious  in  the  case  of  the  first  two  criteria. 
The  choice  of  the  streamline  pattern,  as  the  third  criteria  followed  by  the 
vertical  motion  patterns,  is  based  on  Reed's  (1950)  conclusion  that  the 
horizontal  advection  contributed  2/3  of  the  day-to-day  change  whereas  the 
vertical  motions  contribute  1/3  of  the  daily  change . 

5.  Ozone  Distribution  in  the  Lower  Stratosphere 

Before  discussing  the  individual  ozone  analyses  it  is  worthwhile 
to  compare  the  twelve-day  mean  ozone  mixing-ratio  distribution  obtained  here 
with  distributions  prepared  by  other  investigators. 

The  twelve-day  mean  ozone  mixing  ratio  computed  for  the  period  29  April 
to  10  May  1967  is  presented  in  figure  5.  It  is  seen  that  maximum  values  at 
southern  latitudes  and  high  altitudes  give  way  to  lower  values  at  northern 
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FIG.  5.  Twelve-day  Zonal  Mean  Ozone  Mixing  Ratio  (mg  g-1)  for  the  period 
29  April  to  10  May  1963.  Computed  between  longitudes  60W  and  140W. 


latitudes.  The  latitudinal  gradient  reverses  between  the  30  and  50  mb  surfaces 
with  the  lower  altitudes ,  having  higher  values  at  northern  latitudes.  The 
latitudinal  gradient  reverses  between  the  30  and  50  mb  surfaces  with  the 
lower  altitudes,  having  higher  values  poleward.  The  slope  of  the  isolines 
between  50  and  200  mb  is  markedly  downward.  In  comparison  the  isolines 
between  50  and  30  mb  are  quasi-horizontal  between  35N  and  50N  and  rise 
northward  between  5 ON  and  7 ON. 

A  direct  comparison  of  ozone  mixing  ratio  values  may  be  made  with  the 
values  presented  by  Hering  (1966)  for  the  spring  of  1963-1964.  The  two 
figures  are  similar  with  the  exception  that  slightly  higher  values  are  given 
hereby  using  the  29  April  1963  to  10  May  1963  data.  Although  Hering  uses  the 
same  data  source  as  this  study,  these  are  the  only  available,  directly 
comparable,  mean  ozone  mixing  ratio  values  computed  from  vertical  soundings. 

Bojkov  (1969)  presents  an  annual  vertical  ozone  distribution  as  determined 
by  the  Umkehr  method.  The  units  of  this  distribution  are  micromillibars  and 
not  easily  comparable  in  a  quantitative  sense.  Qualitatively,  the  same 
conclusions  revealed  by  the  twelve-day  mean  values  are  observed  in  Bojkov's 
distribution. 

Taking  into  account  the  information  presented  by  Bojkov  and  Hering,  it  is 
concluded  that  the  twelve-day  mean  cross  section  is  in  agreement  with  their 
results,  and  therefore,  lends  credibility  to  the  individual  analyses  from  which 
the  mean  distribution  was  prepared.  A  representative  sample  of  the  distribution 
of  ozone  mixing  ratio  on  the  200,  100,  50  and  30  mb  pressure  surfaces  is 
presented  in  Appendix  A. 
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The  p<  ‘terns  made  by  the  30  mb  ozone  mixing  ratio  isolines  are  dominated 
by  the  effects  of  the  horizontal  component  of  the  wind.  Ozone-rich  air  is 
associated  with  the  northward  flow  of  the  subtropical  anticyclones  forming 
on  the  southern  edge  of  a  weakening  winter  circumpolar  vortex.  Ozone-poor 
air  is  associated  with  the  southward  flow  on  the  western  side  of  the  large 
scale  troughs  superimposed  on  the  circumpolar  vortex.  The  two  contrasting 
air  masses  meet  in  the  middle  latitudes  and  mix.  The  air,  continuing  to  move 
eastward,  then  diverges  with  relatively  ozone-poor  air  associated  with  the 
southward  flow  on  the  eastern  side  of  the  subtropical  anticyclones.  Corre¬ 
spondingly,  relatively  ozone-rich  air  is  associated  with  the  northward  flow 
on  the  eastern  side  of  the  large  scale  troughs . 

If  long  meridional  trajectories  form  in  connection  with  these  cyclonic  and 
anticyclonic  wind  systems,  unusually  large  values  of  ozone  can  be  observed 
at  northern  latitude  stations  while  lower  values  are  observed  several  hundred 
miles  cast  or  west  of  that  station.  The  converse  holds  true  for  southern 
latitude  stations.  The  movement  of  these  wind  systems  will  affect  the  daily 
observations  at  a  station  in  that  large  day-to-day  fluctuations  of  ozone  can  be 
caused  by  the  associated  wind  shift. 

By  comparing  the  vertical  motion  fields  on  a  pressure  surface  (not  repro¬ 
duced  here)  with  the  ozone  mixing  ratio  distribution,  some  very  interesting 
observations  can  be  made.  The  western  side  of  the  anticyclones  are  associated 
with  rising  vertical  motion,  northerly  trajectories,  and  high  ozone  values 
while  the  western  side  of  troughs  are  associated  with  sinking  motions, 
southerly  trajectories,  and  low  ozone  values.  Therefore,  it  is  qualitatively 
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observed  that  the  Reed-  Normand  effect  is  operating  in  the  opposite  sense 
than  at  lower  altitudes.  The  result  of  this  mechanism  is  in  evidence  in  the 
mean  ozone  mixing  ratio  distribution  of  figure  5  where  the  isolines  between 

•j 

00  and  50  mb  slope  upward  at  northern  latitudes. 

Masses  of  air  with  similar  ozone  characteristics  can  be  identified  and 
tracked  for  periods  of  several  days.  Such  identifiable  air  masses  seem  to  be 
prevalent  in  the  wind  fields  of  the  southern  anticyclones. 

Ozone  mixing  ratio  patterns  at  50  mb  are  similar  to  the  patterns  at  30  mb 
with  ..everal  exceptions.  Streamline  patterns  at  this  level  reflect  ti  e  upper, 
weakening  portions  of  tropospheric  systems.  Weaker  southern  anticyclones 
are  also  present.  The  ozone  mixing  ratio  pattern  shows  ozone-rich  air 
associated  with  the  southward  flow  that  at  30  mb  is  associated  with  ozone- 
poor  air.  Similarly,  ozone-poor  air  is  associated  with  the  northward  flow  that 
at  30  mb  is  associated  with  ozone-rich  air. 

A  qualitative  comparison  of  the  vertical  motion  patterns  with  the  ozone 
mixing  ratio  patterns  show  the  Reed-Normand  effect  to  be  operating  in  the 
manner  originally  proposed  by  Reed  (1950)  and  Normand  (1953).  The  sinking, 
southward  moving  air  into  troughs  are  areas  of  high  ozone  values  while  rising 
northward  moving  air  is  associated  with  low  ozone  values. 

Masses  of  ozone-rich  and  ozone-poor  air  can  be  identified  as  at  30  mb. 
However,  at  both  pressure  levels  it  noted  that  these  identifiable  ozone  masses 
may  be  an  unnatural  consequence  of  the  trajectory  supplementation  of  the 
data  and  a  weak  latitudinal  gradient  of  ozone. 

The  ozone  mixing  ratio  patterns  at  100  and  200  mb  are  similar  to  the  50 
mb  patterns.  Tne  effect  of  the  streamline  pattern  is  more  pronounced. 
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The  vertical  motion  pattern  on  a  borad  scale  reflects  the  Reed-Normand 
effect;  however,  vertical  motion  maxima  with  horizontal  dimensions  on  the 
order  2  or  3  grid  intervals  are  not  clearly  associated  with  the  Reed-Normand 
predicted  high  and  low  ozone  values.  This  is  caused  by  the  density  of  the 
ozone  data,  being  insufficient  to  isolate  these  smaller  scale  features  or  by 
bad  a;  calculations.  Support  for  this  point  of  view  is  gained  from  a  calcu¬ 
lation  of  the  local  change,  horizontal  advection  and  vertical  advection  terms 
of  the  Eulerian  expansion  of  the  total  derivative  of  ozone  mixing  ratio  with 
respect  to  time.  This  calculation  was  made  on  a  point  by  point  basis  for 
the  entire  grid.  The  results  were  inconclusive  since  the  error  term,  computed 
by  assuming  ozone  is  conserved  and  summing  the  three  remaining  terms,  was 
too  large  in  most  cases.  This  indicated  a  smaller  scale  ozone  distribution 
that  was  r.ot  able  to  be  identified  with  the  data  available. 

The  ozone  analyses  prepared  for  this  study  compare  in  a  general  manner 
with  the  distributions  used  by  Breiland  (1964,  1965,  1967,  1968),  Sticksel 
(1965),  Pittock  (1969),  Briggs  and  Roach  (1963),  Danielsen  (1968),  and 
Berggren  and  Labltzke  (1965,  1968)  in  their  respective  studies.  Only  the 
200  mb  analyses  for  3  and  4  May  1963  presented  by  Berggren  and  Labitzke 
(1968)  are  directly  comparable  with  the  analyses  compiled  here.  The  analyses 
are  similar  in  that  the  major  features  are  the  same;  yet  minor  differences 
occur  with  the  smaller  features.  It  is  significant  to  note  that  although  the 
AFCRL  data  was  the  initial  input  in  both  cases,  Berggren  and  Labitzke  used 
isentropic  trajectories  to  supplement  their  data  points  while  this  study  used 
24-hour  isobaric  trajectories  corrected  for  vertical  motion. 
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With  the  data  apparently  consistent  with  the  earlier  theories,  previous 
investigators,  and  consistent  with  respect  to  other  atmospheric  parameters, 
it  was  decided  that  these  analyses  would  provide  an  acceptable  data  base 
for  a  quantitative  calculation  of  ozone  transport  mechanisms. 

E.  THE  OZONE  BALANCE  EQUATION  FOR  A  VOLUME  BOUNDED  BY  TWO 
PRESSURE  SURFACES 

As  a  first  step  in  evaluating  earlier  theories  of  ozone  distribution  and 
movement,  the  calculation  of  the  mean  and  eddy  fluxe  of  ozone  on  a  synoptic 
scale  is  desirable.  Considering  the  minimum  amoun  oi  ozonesonde 
observations  at  a  given  time,  the  uncertainty  in  the  synoptic  distribution  of 
ozone  on  a  pressure  surface,  and  the  uncertainty  in  the  performance  of  the 
ozonesonde,  an  averaging  scheme  over  many  data  points  is  preferred  as 
compared  to  a  point  by  point  evaluation.  To  achieve  this,  the  continuity 
equation  for  ozone  is  written  in  pressure  coordinates  as 


'ff'-V’x- “0+  (p-R)+t;  (Kpf£>  +’-Vx  (3) 

By  putting  the  advection  terms  in  equation  (3)  into  flux  form  and  taking 
the  area  average,  equation  (3)  becomes 

-»TXW2-f^+  (p-R>+fj  kpi?-1-  •  (4) 

Now  by  applying  the  divergence  theorem ,  equation  (4)  may  be  rewritten 


as 


a*  =I 

a  t  a 


j>  x 


c  dl 
n 


_ 

dp 


+  (P  -  R) 


+  f-  <K  !*) 
3  p  p  3  p 


K  I* 

H  3n 


dl. 


(5) 


where  c  ,  the  component  of  the  wind  normal  to  the  boundary,  is  defined  as 
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positive  toward  the  center  of  the  volume. 


Along  a  line  let  x  =  X'*  where  x  is  the  line  average  and  x'  is  the 
point  deviation  from  the  line  average.  Similarly,  for  averaging  operations, 
equation  (5)  becomes 

+  7  Tc-1-  .  l22. 

3 1  A  A  *  n  3p  dp 

<- - -  K  _  (6) 

+  (P  -  R)  +  (k  !*  )  -  —  f*  A 1  , 

3p  \  p  3p  /  A  3n 

where  it  is  assumed  that  K-,  is  constant  on  a  pressure  surface  and  noted 

H 

that  j>  dl  -  A  1 . 

Now  average  equation  (6)  between  two  pressure  surfaces,  p^  and 
P^,  to  obtain  the  ozone  balance  equation 


(LC)  (MFLA) 


(EFLA) 


3y  A1  -  -  A  1  — - r 

3t  A  *  n  A  *  n 


(MFLO)  (MFT) 

<p_X  |  +  co  y  | 

AP  'pl  AP  >pu 


(EFLO) 


(EFT) 


Ap  'PT  Ap  'p 


(PR)  (SUBLO) 

/VVV^  ^ 

+  (P^R)  +  | 


Ap  3p  PT 


(SUBT) 


(SUB  LA) 


_  _£  ix  |  _  M  ( K  \ 

Ap  ap  'p  A  \  H  3n  /  ' 

where  the  vertical  average  has  been  applied  to  evaluate  the  terms  in  a  form 
for  numerical  computations. 

The  terms  of  the  equation  are: 


(LC)  the  time  rate  of  change  of  average  ozone  mining  ratio 
over  the  volume; 
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(MFLA)  the  mean  flux  of  ozone  into  the  region  through  the  lateral 
boundaries  (the  vertical  side  of  the  volume); 

(EFLA)  The  eddy  flux  of  ozone  through  the  lateral  boundaries; 

(MFLO)  the  mean  flux  of  ozone  into  the  region  across  the 
bottom  boundary; 

(MFT)  the  mean  flux  of  ozone  into  the  region  across  the  top 
boundary; 

(EFLO)  the  eddy  flux  of  ozone  into  the  region  across  the  bottom 
boundary; 

(EFT)  the  eddy  flux  of  ozone  into  the  region  across  the  top 
boundary; 

(PR)  the  average  difference  between  production  and  recombin¬ 
ation  of  ozone  for  the  whole  volume; 

(SUBLO)  the  flux  of  ozone  into  the  region  across  the  upper  boundary, 
due  to  sub-grid-scale  diffusion; 

(SUBLA)  the  flux  of  ozone  into  the  region  through  the  sides  due  to 
sub-grid-scale  diffusion. 

All  the  terms  have  units  of  micrograms  per  gram  per  day.  This  unit  is 
convenient  in  view  of  the  time  scale  of  the  observations.  In  addition,  the 
value  of  the  ozone  mixing  ratio  at  a  given  height  appears  as  a  variable  on  the 
ozonagram. 

Computer  programs  were  written  to  evaluate  all  terms  on  a  daily  basis  for 
the  30/50,  50/100,  and  100/200  mb  volumes  with  the  exception  of  the  produc¬ 
tion  and  recombination  term  and  the  sub-grid  scale  diffusion  terms . 

Brewer  and  Wilson  (1968)  have  evaluated  the  meridional  net  ozone 
production  rate  for  the  winter  and  summer  season  from  the  equator  to  sixty 
degrees  latitude  and  between  200  and  30  mb.  Since  the  evaluation  of  the  PR 
term  would  require  a  complete  and  separate  treatment.  Brewer  and  Wilson's 
results  are  accepted  for  use  in  this  calculation.  To  compute  an  estimated 
value  and  the  net  ozone  production  rate,  a  rough  interpolation  between  seasons 


III—  3  7 


and  pressure  surfaces,  followed  by  an  estimated  latitudinal  average,  was 
necessary.  The  appropriate  unit  conversions  were  made.  These  values 
will  be  presented  with  the  results  of  the  ozone  balance  equation  In  a  later 


section. 

The  sub-grid-scale  diffusion  terms  are  expected  to  be  small  when 
compared  to  the  magnitude  of  the  other  terms.  Lacking  a  reliable  estimate 
of  the  eddy-diffusion  coefficients  of  ozone  mixing  ratio,  these  terms  were 
not  calculated. 

The  initial  attempt  to  evaluate  the  ozone  balance  equation  by  using  the 
observed  wind  to  directly  evaluate  each  term  resulted  in  unusually  large 
errors  that  could  be  traced  to  the  fluctuations  of  the  mean  terms.  To  improve 
upon  these  results,  the  mass  continuity  equation  was  solved  for  the  mean 
normal  component  of  the  wind  at  the  side  boundaries  and  in  terms  of  the 
computed  to  fields .  The  evaluation  of  the  mean  normal  component  of  the 
wind  can  be  obtained  by  applying  the  divergence  theorem  to  the  area-averaged 


continuity  equation  and  solving  the  following  form, 

a  -  w  ) 


C  = 


PL  "  ?u 


n  A1  (pL  -  pu) 


(8) 


Although  this  scheme  assured  the  conservation  of  mass,  it  increases  the 
dependence  of  the  calculations  on  the  computed  vertical  motions  as  compared 
to  directly  observed  winds. 
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III.  EVALUATION  OF  THE  TERMS 
OF  THE  OZONE  BALANCE  EQUATION 


Without  further  modification  the  terms  of  the  ozone  balance  equation 
were  computed.  The  daily  calculations  tended  to  fluctuate  rapidly  and  may 
be  unrealistic.  Considering  the  magnitude  of  the  error  term  in  comparison 
with  the  magnitude  of  the  other  terms,  it  was  decided  to  present  the  ten-day 
mean  values  of  the  terms.  These  values  are  given  in  tabular  format  in 
Table  II.  In  the  computed  format,  the  individual  mean  and  eddy  terms  cannot 
be  compared  directly.  To  effect  a  comparison  the^mean  terms  must  be  summed, 
thereby  giving  the  net  effect  of  the  mean  terms  on  the  volume. 

To  obtain  a  form  in  which  the  mean  horizontal  and  vertical  flux  may  be 
directly  compared  to  the  eddy  flux  terms,  the  sum  of  the  mean  terms,  as 
presented  in  Table  II,  may  be  computed  in  another  manner.  The  area- 
averaged  continuity  equation  can  be  written  as  a  representation  of  the  net 
effect  of  the  mean  terms  at  the  middle  pressure  between  two  pressure 
surfaces.  By  applying  a  vector  identity  and  the  divergence  theorem,  the 
first  and  third  terms  on  the  right  hand  side  of  equation  (6)  can  be  written  as 

SUM  =  ~  x  ~  -  i  -  J  i  e  6  1  .  (9) 

n  *  A  op  A  J  n  '  ' 

By  combining  like  terms  and  putting  into  computational  format,  equation 

(9)  becomes 


(MLAD)  (MVAD) 


(10) 
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Key:  See  Section  II-E 


SUM  =  MFLA  +  MFLO  +  MFT 
AS  COMPUTED: 

Volume  (LC)  (MFLA)  (EFLA)  (MFLO)  (MFT)  (EFLO)  (EFT  (PR)  (ERROR) 


30/50  mb 

-173  =  +  184  +  54 

-688 

+696 

-  67 

-50 

-125  +177 

50/100  mb 

-517  =  -1102  -  46 

-921 

+275 

-103 

+27 

-  44  +806 

100/200  mb 

-273  =•  -  229  -321 

-  47 

+461 

-155 

+51 

-16+17 

COMPARISON  FORMAT: 


Volume  (LC) 

(SUM) 

(EFLA) 

(EFLO) 

(EFT) 

(FR) 

(ERROR) 

30/50  mb  -173 

=  +192 

+  54 

-  67 

-50 

-125 

+  177 

50/100  mb  -517 

=  +456 

-  46 

-103 

+27 

-  44 

+806 

100/200  mb  -273 

-  +185 

-321 

-155 

+51 

-  16 

+  17 

TABLE  II.  The  Ten-day  Mean  Ozone  Balance  for  Indicated 
Volumes  .  Symbols  for  each  term  given  in 
Section  II— E .  Units  are  micro', rams  per  gram 
4 

per  day  x  10  . 


III-40 


The  two  terms  are  defined  as: 


(MLAD)  the  mean  advection  through  the  later  boundary; 

(MVAD)  the  advection  into  the  volume  due  to  the  mean 

vertical  motion . 

The  two  terms  of  equation  (10)  were  evaluated,  using  the  same  input 
values  as  for  the  evaluation  of  equation  (7).  The  results  of  the  evalua¬ 
tion  of  equation  (10)  are  presented  in  pictorial  form  in  figure  6  with  the 
local  change  and  eddy  terms  as  evaluated  by  equation  (7).  The  net 
production  term  is  the  same  in  both  cases.  The  error  terms  in  each  evalua¬ 
tion  are  different  only  to  a  small  degree.  It  is  therefore  concluded  that 
the  computational  scheme,  used  to  arrive  at  the  value  of  the  mean  terms,  is 
a  reasonable  evaluation  of  the  input  parameters. 

With  the  mean  and  eddy  ozone  flux  terms  in  a  format  for  comparison  and 
some  degree  of  confidence  in  the  method  of  evaluation  of  the  mean  terms  , 
it  is  appropriate  to  discuss  the  individual  terms  of  the  ozone  balance  as 
presented  in  figure  6. 

A.  THE  LOCAL  CHANGE  TERM 

Computation  of  this  term  was  made  directly  from  tne  distribution  of 
ozone  on  a  pressure  surface,  and  therefore ,  is  not  affected  by  other  computed 
parameters.  A  rough  check  on  the  value  of  this  term  was  made  by  comparing 
the  computed  result,  multiplied  by  100,  with  a  100-day  difference  between 
area  averages  estimated  from  Godson's  (1960)  mean  distribution  of  total  ozone. 
This  period  falls  after  the  March  maximum  of  total  ozone  shown  by  Godson  and 
the  sign  of  the  mean  local  change  term  Is  negative.  With  these  observations 
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FIG.  6.  Ten-day  mean  ozone  balance  in 
advective  format.  Symbols  for  each  term 
Riven  in  Section  II-E  and  Section  IIT.  , 
Units  are  microRrams  per  Rram  per  dav  x  10 


in  mind,  the  100-day  period  between  May  and  August  was  chosen  for 
comparison . 

After  making  the  appropriate  unit  conversions,  the  mean  local  change 
term  accounted  for  85  percent  of  the  change  in  Godson's  distribution.  Thus, 
the  rate  of  decrease  of  ozone  within  the  volume  for  this  period  appears 
representative  of  the  seasonal  trend. 

It  is  significant  to  note  that  the  maximum  of  total  ozone  in  the  atmosphere 
is  observed  in  March  and  a  net  hemispheric  decrease  is  observed  after  that 
month  (Godson,  1960).  The  dynamic  mechanisms  or  conditions  that  created 
the  maximum  conditions,  particular’y  at  the  northern  latitudes,  must  have 
ceased  operating.  One  may  speculate  that  the  mechanism  to  accomplish  this 
winter-spring  northern  maximum  is  the  long  meridional  trajectories  set  up 
by  the  breakdown  of  the  polar  night  vortex  associated  with  the  mid-winter 
"sudden  warmings."  Further,  the  dynamic  mechanisms  or  conditions  now 
operating  must  be  serving  to  deplete  the  atmosphere  of  ozone. 

B.  THE  MEAN  TERMS 

The  u-;an  advection  through  the  lateral  boundaries  can  only  be  evaluated 
as  to  direction  into  or  out  of  the  volume.  In  view  of  this  limitation,  the 
only  point  to  be  made  is  a  mention  that  the  term  in  each  volume  is  contributing 
a  small  net  loss  to  the  volume  and  the  mean  terms  (MLAD)  are  smaller  than 
their  respective  eddy  terms  (EFLA)  (see  fig.  6). 

On  the  other  hand,  the  mean  advection  in  the  vertical  shows  a  consistent 
downward  motion  through  the  three  volumes  of  interest.  In  the  two  upper 
volumes,  the  mean  vertical  advection  is  the  dominant  term. 
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C.  THE  EDDY  TERMS 


In  the  30/50  mb  volume  the  lateral  eddy  flux  term  (EFLA)  Is  an  input  to 
the  volume  and  agrees  with  the  earlier  qualitative  observation  that  the 
anticyclones  on  the  southern  edge  of  the  circumpolar  vortex  appeared  to  be 
transferring  ozone  northward.  The  magnitude  of  the  lateral  eddy  flux  term 
(EFLA)  is  small  relative  to  the  observed  change  within  the  volume. 

The  lateral  eddy  flux  terms  at  50/100  and  100/200  mb  serve  to  deplete 
the  volumes  of  ozone.  The  value  of  the  term  at  50/100  mb  is  small  relative 
to  the  mean  vertical  advection  and  the  vertical  eddy  flux  across  the  100  mb 
surface.  In  the  100/200  mb  volume  the  lateral  eddy  flux  term  is  the  largest 
term  in  the  ozone  balance.  An  interesting  and  unexplained  question  is  the 
ultimate  destination  of  the  ozone  advected  out  of  the  volume  of  computations 
by  the  lateral  eddy  flux. 

The  vertical  eddy  flux  terms  show  a  consistent  downward  transport  from 
50  to  200  mb  and  are  a  measure  of  the  vertical  motion  portion  of  the  Reed- 
Normand  effect.  At  these  levels  the  direction  is  consistent  with  that  pro¬ 
posed  by  Reed  and  Normand.  At  the  30  mb  surface  the  direction  of  the 
vertical  eddy  flux  is  upwards.  This  can  be  traced  to  the  northward,  rising 
ozone-rich  air,  associated  with  the  western  side  of  the  anticyclones  at  this 
level.  This  is  consistent  with  the  e^.iier  section  in  which  the  Reed-Normand 
effect  was  observed  to  be  operating  in  the  opposite  sense  as  observed  in 
lower  altitudes  . 

By  comparing  the  vertical  eddy  flux  at  the  upper  and  lower  surfaces 
at  each  volume,  it  can  be  seen  that  these  eddies  contribute  a  net  loss  of 
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ozone  to  each  volume.  Another  interesting,  unanswered  question  is  posed 
by  this  observation.  If  a  given  amount  of  ozone  is  going  downward  into  a 
trough,  what  physical  process  or  dynamic  mechanisms  are  operating  to 
prevent  an  equal  amount  of  ozone  from  proceeding  back  upwards  to  the  ridge? 
Tie  answer  probably  lies  with  a  smaller  scale  turbulence  phenomenon  not  yet 
identified  much  as  the  energy  cascade  in  the  lowest  level  of  the  atmosphere. 

The  vertical  eddy  flux  terms  are  the  most  significant  of  all  the  terms 
computed  in  the  ozone  balance.  This  conclusion  is  attributed  to  the  fact 
that  the  vertical  eddy  flux  computed  for  this  particular  area's  pressure  surfaces 
may  be  representative  of  the  whole  global  pressure  surface  in  question.  This 
cannot  be  said  of  the  lateral  eddy  flux  terms.  As  is  to  be  described  in  a 
later  section,  the  magnitude  of  the  mean  terns  is  subject  to  question. 

D.  THE  PRODUCTION  AND  RECOMBINATION  TERM 

In  comparison  with  the  magnitudes  of  the  other  terms,  the  production 
and  recombination  term,  as  estimated  from  Brewer  and  Wilson’s  (1968)  values, 
does  not  have  a  significant  effect  except  in  the  30/50  mb  volume.  Here, 
the  magnitude  is  about  E/3  of  the  local  rate  of  change.  In  future  investigations 
this  must  be  taken  into  consideration. 

E.  THE  ERROR  TERM 

The  error  term  in  Table  II  and  figure  6  is  the  value  that  must  be  added 
to  the  local  change  in  order  to  equal  the  sum  of  the  mean,  eddy  and  net 
pioduction  terms.  In  numerical  computations  of  this  sort,  errors  are  expected 
because  of  data  uncertainties,  since  many  approximations  are  made  to  facil¬ 
itate  the  evaluations  of  the  individual  terms  and  their  components . 
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The  mean  errors  in  each  volume  are  of  the  order  of  magnitude  of  the 
largest  term.  The  mean  errors  at  30/50  and  50/100  mb  are  representative 
of  the  daily  errors  at  these  levels  from  which  the  mean  error  was  derived. 

This  is  not  true  at  the  100/200  mb  volume  where  the  daily  errors  average 
one  order  of  magnitude  larger  than  the  ten-day  mean  error.  Thus,  it  is 
noted  that  the  ten-day  mean  error  is  unrepresentative  of  the  error  in  the 
100/200  mb  volume  and,  further,  that  the  calculations  here  suffer  from  the 
same  problems  as  the  upper  two  volumes. 

Considering  that  the  error  in  all  three  volumes  indicates  too  much 
ozone  on  the  right  hand  side  of  the  balance,  the  failure  to  compute  the  sub¬ 
grid  scale  diffusion  terms  has  to  be  considered.  It  is  expected  that  these 
terms  will  evacuate  ozone  from  the  volume  with  a  magnitude  considerably 
smaller  than  the  magnitude  of  the  error.  Therefore  it  is  improbable  that  the 
failure  to  calculate  these  terms  contributes  much  to  the  error. 

The  first  serious  possible  source  of  error  is  that  the  ozone  mixing  ratio 
distributions  on  the  pressure  surfaces  are  not  accurate  enough  for  the  calcu¬ 
lations.  Only  by  having  a  more  dense  data  network  will  this  ever  be  known 
for  sure.  The  probability,  that  this  is  the  major  source  of  error,  is  diminished, 
considering  the  consistency  in  daily  analyses  and  the  consistency  of  the 
distributions  with  the  theories  and  distributions  advanced  by  earlier  inves¬ 
tigators  . 

The  second  and  most  probable  source  of  error  is  in  the  evaluation  of  the 
vertical  motions.  As  was  noted  earlier,  the  evaluation  of  the  ozone  balance 
equation  is  heavily  dependent  upon  vertical  motions.  While  a  consistent 


error  in  the  vertical  motion?  wii;  not  affect  aif  eddy  terms  greatly,  it  will 
have  considerable  effect  on  the  mean  terms.  Since  the  continuity  equation 
was  used  to  evaluate  the  mean  flux  through  the  lateral  boundaries,  this 
term  is  affected  also.  This  line  of  reasoning  is  consistent  with  the  fact 
that  the  mean  terms  in  most  cases  are  large  compared  to  the  rest  of  the  terms 
as  can  be  seen  in  Table  li  and  figure  6.  It  can  be  reasoned  that  a  consis¬ 
tent  error  in  the  vertical  motion  at  a  particular  pressure  level  could  affect 
each  of  the  mean  terms  enough  to  produce  a  significant  error. 

The  size  of  the  error  in  the  vertical  motion  modes  does  not  have  to  be  very 

large  to  seriously  affect  the  calculations.  To  see  this,  consider  the  vertical 

motion  scheme  used  in  the  computations  and  presented  as  equation  (1).  The 

diabatic  heating  term  is  the  term  most  subject  to  question  as  its  values 

determined  by  any  method  are  highly  uncertain.  This  is  attributed  to  the 

lack  of  an  accurate  method  to  check  the  results  of  the  determination  of  the 

diabatic  heating.  If  the  area  average  of  equation  (1)  is  taken,  and  it  is 

assumed  that  the  correction  to  the  vertical  motion  (a*  )  is  a  function  only 

corr 

of  the  correction  to  the  diabatic  heating  term  (corr),  then  it  can  be  written 
that  approximately, 

icorr‘^E i_  .  in) 

o/cp 

A  representative  value  of  area  averaged  diabatic  heating  at  75  mb  as 
estimated  from  Kennedy's  (1965)  results  can  be  considered  -.35  degrees  per 
day.  If  this  were  incorrect  by  an  arbitrarily  selected  +  .2  degrees  per  day, 
then  the  evaluation  of  equation  (11),  using  representative;  values  for  evalu¬ 
ation  of  the  specific  volume  would  produce 
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U)  =  +  .24  nib  per  day,  or  about  a  25  percent  error  in  the  mean 
con 

vertical  motions . 

The  mean  vertical  advection  of  ozone  across  a  pressure  surface  is 
represented  by  Zi  in  equation  (10).  Evaluating  the  effect  of  ^corr 
on  this  term  using  representative  values  of  x  produces  uncertainty  in 
this  term  of  approximately  210  x  10  *  micrograms  per  gram  per  day.  It 
is  easily  seen  from  equations  (7)  and  (8)  how  this  would  alfect  the  mean  flux 
through  the  lateral  boundaries  as  presented  in  Table  II.  The  value  of  the 
mean  advection  through  the  lateral  boundary  is  small  (Fig.  6).  Therefore, 
it  is  reasoned  that  in  this  form  the  mean  lateral  advection  is  not  a  major 
source  of  error.  Similar  arguments  can  be  made  for  the  30/50  and  100/200 
mb  mean  terms  , 

It  can  be  concluded  that  the  uncertainty  in  the  evaluation  of  mean 
stratospheric  vertical  motions  is  large  enough  to  result  in  the  observed 
magnitude  of  the  error  to  the  ozone  balance  equation  through  evaluation  of  the 
mean  terms.  Implied  from  the  excess  of  ozone  represented  by  the  error  term 
is  that  the  net  input  to  the  particular  volume  by  the  mean  terms,  particularly 
the  mean  vertical  flux  (or  advection)  terms,  is  too  large  although  the  sense 
of  the  individual  mean  terms  is  probably  correct. 
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IV.  ADDITIONAL  CALCULATIONS 


An  attempt  was  made  to  support  the  results  of  the  qualitative  evaluation 
of  the  distributions  of  ozone  on  a  pressure  surface  and  the  terms  of  the  ozone 
balance  by  calculating  the  correlation  coefficients  of  the  meridional  wind 
and  ozone  mixing  ratio  and  the  correlation  coefficients  of  the  vertical  motion 
and  the  ozone  mixing  ratio.  The  calculations  are  presented  as  figures  7  and 
8.  For  the  number  ot  data  points  involved,  a  correlation  coefficient  greater 
than  .2  is  statistically  significant  at  the  95  percent  level  of  confidence. 

The  correlation  coefficient  for  two  variables  is  defined  as  the  covariance 
divided  by  the  variance.  It  follows  then  that  the  sign  of  the  correlation 
coefficient  is  also  the  sign  of  the  covariance.  This  relationship  allows 
comparison  with  the  eddy  terms  of  the  ozone  balance. 

Looking  at  figure  7,  it  can  be  seen  that  there  are  non-significant  values 
of  the  lateral  transport  computed  by  the  ozone  balance.  The  ozone  distribu¬ 
tions  observed  on  the  pressure  surfaces  are  evaluated  as  having  been  spread 
out  by  the  weak  wind  fields  at  these  levels  and  accomplishing  little  transport. 
The  positive  area  at  30  mb  and  70N  latitude  is  not  conclusive,  since  this 
area  was  determined  by  one  data  point.  It  is  consistent,  however,  with  the 
assumption  that  ozone-poor  air  is  advected  southward  by  the  horizontal  eddy 
flow  at  this  altitude  . 

At  the  100  mb  level  the  negative  values  to  the  north  are  associated  with 
high  ozone  values  being  moved  southward  by  the  horizontal  eddy  flow.  The 
positive  values  to  the  south  are  an  indication  of  a  counter-gradient  transport 
(see  fig.  5)  of  high  ozone  values  northward  and  low  ozone  values  southward 
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by  the  horizontal  eddy  flow.  To  evaluate  the  cause  of  the  counter-gradient 
transport,  the  vertical  motion  patterns  (not  presented  here)  were  compared  with 
the  streamline  patterns,  it  was  observed  that  in  many  cases  the  maximum 
sinking  areas  connected  with  ozone -rich  air  were  slightly  (45  degrees  in 
phase  angle)  to  the  east  of  the  base  of  the  trough.  This  offset  of  the  vertical 
motion  pattern  from  the  trough  line  has  been  discussed  by  Ohring  and  Muench 
(1960).  The  sinking  consequently  offset  the  entire  ozone  mixing  ratio  pattern 
to  the  east  and  set  the  conditions  for  counter-gradient  flow  by  the  horizontal 
eddies . 

The  results  at  200  mb  were  expected  to  be  more  conclusive  than  shown  in 
figure  7.  The  values  reflect  the  stronger  relationships  at  100  mb,  however, 
they  are  not  at  the  95  percent  level  of  significance .  Without  further  calcula¬ 
tions ,  a  definite  conclusion  cannot  be  advanced. 

The  results  of  figure  8  strongly  support  earlier  contentions.  The  weak 
positive  values  are  associated  with  upward  eddy  transport  of  high  ozone  values 
by  the  anticyclonic  activity  at  3n  and  50  mbs.  The  high  negative  values  at 
the  southern  latitudes  indicate  the  exoected  downward  transport  of  high  ozone 
values.  These  latitudes  coincide  with  the  base  of  the  troughs  on  the  circum¬ 
polar  vortex  and  hence  general  areas  of  maximum  sinking  motions.  At  100  mb 
the  secondary  maximum  o,  significant  negative  values  at  the  northern  latitudes 
are  evaluated  as  downward  transport  of  high  ozone  values  by  the  weaker  semi¬ 
permanent  trough  of  the  circumpolar  vortex.  Based  on  the  values  at  200  mb, 
it  is  further  reasoned  that  the  vertical  eddy  flux  of  ozone  across  the  tropopause 
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is  talcing  place  in  connection  with  the  maximum  sinking  areas  near  the  base 
of  the  transient  troughs  on  the  circumpolar  vortex. 

In  the  recent  past  several  investigators  (Newell  1961,  1962,  1963,  1964; 
Mahlman,  1966,  1967;  and  Molla  and  Loisel,  1962)  have  used  the  latitudinal 
covariance  of  the  meridional  and  vertical  wind  components  (v'w')  as  an 
indication  of  stratospheric  eddy  transport  processes.  As  a  further  check 
on  the  validity  of  the  terms  of  the  ozone  balance  equation,  it  was  decided  to 
compare  (v'w1),  by  using  correlation  coefficients,  with  the  information  prepared 
for  this  study.  This  distribution  is  presented  in  figure  9. 

The  negative  values  of  the  correlation  coefficients  are  associated  with  a 
northward  and  downward  transport  while  the  positive  areas  indicate  a  south¬ 
ward  and  downward  transport.  The  information  from  fig.  9  is  in  agreement  with 
the  sense  of  the  vertical  eddy  flux  terms  of  the  ozone  balance,  and  the  ozone 
distribution  on  a  pressure  surface  except  at  the  30  mb  surface.  At  this  level 
the  mean  ozone  mixing  ratio,  streamline  and  vertical  motion  patterns  show  a 
northward  and  upward  flow  of  high  ozone  values.  This  was  confirmed  by  a 
vertical  eddy  flux  out  of  the  30/50  mb  volume  across  the  30  mb  surface  and 
figure  8.  The  vertical  eddy  flux  is  a  direct  calculation  ox  transport,  as 
compared  to  the  correlation  coefficients  which  are  only  an  indication  of  trans¬ 
port.  This  implies  that  the  correlation  coefficient  must  be  on  the  order  of 
magnitude  of  .5  or  greater  to  give  confidence  in  the  conclusion  that  this  term 
is  in  lact  an  indicator  of  transport. 


longitude  60W  and  140W. 


V.  SUMMARY  AND  CQKCl.USTONS 


Even  tnouqh  the  magnitude  of  tlie  mean  terms  is  subject  to  question,  the 
distribution  of  ozone  on  a  pressure  surface,  the  ozone  balance,  and  the 
supporting  calculations  allow  several  general  conclusions  to  be  drawn. 

1.  It  is  possible  to  construct  the  distribution  of  ozone  on  a  pressure 
surface  from  the  Regener  ozonesonde  uata,  and  then  perform  numerical  calcu¬ 
lations  that  provide  reasonable  results. 

2.  During  this  time  period,  the  atmosphere  is  past  its  yearly  maximum 
content  of  ozone  at  northern  latitudes.  There  is  no  evidence  to  support  a  large 
northward  movement  of  ozone  at  this  time . 

3.  During  the  ten-day  period  in  question,  ozone  content  of  the  volume 
is  decreasing.  To  accomplish  this,  the  mean  vertical  motions,  supported  by 
the  vertical  eddy  transport  at  50,  100,  and  200  mb,  are  moving  ozone  downward 
to  the  tropopause  level  where  the  vertical  eddy  transport,  associated  with  the 
transient  eddies  on  the  circumpolar  vortex,  then  move  the  ozone  into  the 
troposphere  for  eventual  destruction. 

4.  The  vertical  eddy  transport  is  a  measure  of  the  vertical  part  of  the 
Reed-Normand  effect.  This  mechanism  is  operating  the  sense  as  proposed 
between  50  and  200  mb.  At  30  mb  the  Reed-Normand  effect  is  operating  in  the 
opposite  sense  with  ozone-rich  au  being  advected  upward. 
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APPENDIX  A 


STREAMLINE  AND  OZONE  MIXING  RATIO  ANALYSIS  FOR 

3  MAY  1967 

This  appendix  presents  a  representative  sample  of  the  distributions  of 
ozone  mixing  ratio  on  a  pressure  surface  in  the  lower  stratosphere.  To 
understand  the  effect  of  the  wind  field  upon  the  ozone  distribution,  the 
appropriate  streamline  pattern  is  presented. 

The  choice  of  data  was  made  to  facilitate  comparison  with  the  patterns 
published  by  Berggren  and  Labitzke  (1968). 
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in  the  Atmosphere 
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I.  Introduction 


The  investigation  of  the  mechanisms  in  the  atmosphere,  which  lead  to 
systematic  large-scale  transport  of  trace  substances,  is  still  a  relatively 
new  research  endeavor.  Ideally,  the  most  successful  experimental  setup 
for  such  investigations  would  involve  the  joint  use  of  meteorological  and 
trace  substance  data.  Normally,  both  types  of  data  should  possess  similar 
spatial  and  temporal  resolutions.  Up  to  the  present  time,  this  has  been  nearly 
impossible  because  of  the  tremendous  expense  and  difficulty  involved  in 
obtaining  numerous  simultaneous  measurements  of  any  trace  substance  in  the 
free  atmosphere  at  the  desirable  elevations.  As  a  result  of  these  difficulties, 
previous  investigators  have  employed  a  variety  of  approaches  for  investi¬ 
gating  problems  of  trace  substance  transport. 

One  of  the  goals  of  the  research  conducted  under  this  contract  for  the  past 
two  years  has  been  to  explore  possible  new  techniques  for  investigating 
transport  phenomena  in  the  atmosphere.  Resulting  from  this,  a  number  of 
different  approaches  have  been  attempted.  Some  of  these  are  "traditional" 
in  the  sense  that  other  investigators  have  employed  them  in  the  past.  Others 
are  "new"  in  that  they  had  not  yet  been  applied  to  the  investigation  of  larae- 
scale  trace  transport. 

Because  a  number  of  techniques  have  been  attempted  here,  it  seems 
appropriate  to  the  author  that  the  final  project  summary  be  directed  toward 
an  overview  of  the  utility  and  applicability  of  these  various  approaches. 
Consequently,  the  usual  project  "review"  of  pertinent  results  will  not  be 
attempted  here,  since  those  results  are  already  published  under  individual 
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titles,  here  and  elsewhere. 

The  app'oach  in  this  paper  will  be  to  list  each  study  performed  during 
the  contract  period,  and  then  comment  on  the  method  employed  and  results 
obtained  on  an  individual  basis.  The  publications  under  this  contract  listed 
in  order  of  subsequent  referral  in  the  text,  are  as  follows: 

A.  Mahlman,  J.  D. ,  1969:  Heat  Balance  and  mean  meridional  circulations 
in  the  polar  stratosphere  during  the  sudden  warming  of  January  1958. 
Monthly  Weather  Review,  97,  534-540. 

B.  Mahlman,  J.  D. ,  1969:  Energetics  of  a  “minor"  breakdown  of  the 
stratospheric  polar  night  vortex.  Journal  of  the  Atmospheric  Sciences, 

26,  1306-1317.  " “  ~  ‘  '  '  . . 

(j .  Mahlman,  J.  D.,  1970:  Eddy  transfer  processes  in  the  stratosphere 
during  major  and  "minor"  breakdowns  of  the  polar  night  vortex. 

Journal  of  Geophysical  Research,  75.  1701-1705. 

D.  Mahlman,  J.  D.,  1970:  On  the  maintenance  of  the  polar  front  jet  stream. 
(Enclosed  as  part  of  this  report.) 

E.  Mahlman,  J.  D.,  1969:  Diagnostic  analysis  of  the  transport  mechanisms 
leading  to  surface  fallout  increases  during  November  1966.  Naval 
Postgraduate  School,  First  Progress  Report  to  the  Atomic  Energy 
Commission  (NPS-5 1 MZ9081B) ,  5-6. 

F.  Mahlman,  J.  D.,  1969:  Long-term  dependence  of  surface  fallout 
fluctuations  of  tropopause-level  cyclogenesis  .  Archiv  fur  Meteorologle , 
Geophyslk  und  Bioklimatologles ,  A18,  299-311  . 

G.  Stender,  R.  H.,  1970:  Numerical  calculations  of  stratospheric  ozone 
transport  (enclosed  as  part  of  this  report). 

H.  Willoughby,  H.  E,,  1970:  A  numerical  simulation  of  the  advective  and 
diffusive  transfer  of  trace  substances  in  the  stratosphere  (enclosed  as 
part  of  this  report). 

II.  Discussion 

Articles  A,  B,  C,  and  D  in  the  preceeding  list  fall  under  the  category 
of  what  shall  be  termed  "dynamical  process"  studies.  These  studies  are  only 
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indirectly  related  to  the  large-scale  trace -substance  transport  problem  in 
that  they  specifically  attempt  to  understand  the  structure  and  dynamics  of 
various  circulation  phenomena.  However,  once  such  circulation  characteristics 
have  been  properly  delineated,  it  becomes  possible  to  construct  inferences  on 
their  transport  properties  as  well.  This  approach  has  been  employed  pre¬ 
viously  by  a  number  of  investigators  and  has  provided  valuable  insights  into 
a  great  variety  of  atmospheric  processes.  Indirectly,  these  studies  have  laid 
a  framework  for  viewing  the  trace-substance  transport  problem  by  helping  to 
isolate  the  relative  importance  of  the  various  features  of  the  atmospheric 
general  circulation.  Conceptually,  it  has  been  found  useful  to  separate 
these  features  into  distinct  categories  such  as:  the  mean  meridional  circu¬ 
lation,  "eddies"  or  waves  of  various  scales,  and  sub-meteorological-grid- 
scale  "eddies"  or  "diffusion."  Such  a  separation  is,  of  course,  not  unique 
nor  are  the  motions  on  various  scales  independent  of  one  another.  An 
example  of  an  alternative  framework  for  viewing  atmospheric  processes  is 
given  in  Article  D  of  tills  discussion. 

Articles  A,  B,  and  C  are  oriented  toward  attempting  to  understand  specific 
circulation  phenomena  in  the  wintertime  lower  stratosphere.  These  studies 
have  concentrated  on  shorter  time  scales,  and  have  been  intentionally  designed 
to  complement  the  studies  performed  in  the  Planetary  Circulations  Project  at 
the  Massachusetts  Institute  of  Technology,  which  have  been  concerned  with 
circulations  on  the  climatic  scale. 

la  both  Articles  A  and  B,  one  of  the  major  goals  was  to  deduce  the  mean 
meridional  circulation  in  order  that  both  the  dynamics  and  transport  properties, 
associated  with  specific  circulation  phenomena,  could  be  more  fully  under¬ 
stood.  This  was  accomplished  by  employing  a  heat  balance  approach.  Such 
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a  method  proves  to  give  excellent  results  as  long  as  the  temperature  changes 
induced  by  the  mean  meridional  circulation  are  significantly  larger  than  the 
present  uncertainty  in  the  specification  of  diabatic  heating  rates. 

A  major  breakdown  of  the  stratospheric  polar  night  vortex  was  the  object 
of  study  in  Article  A,  while  Article  B  was  concerned  with  a  so-called  "minor 
breakdown"  event.  In  the  "minor  breakdown”  work,  considerable  effort  was 
employed  to  explain  the  kinetic  energy  variations  occurring  in  association 
with  the  breakdown.  Even  though  quite  consistent  mean  meridional  circulations 
were  obtained,  it  still  was  apparent  that  the  kinetic  energy  balance  is  sensi¬ 
tive  to  small  errors  in  the  mean  cell. 

Article  C  was  a  combination  of  the  research,  lending  to  Articles  A  and  B 
to  obtain  a  comparison  of  the  relative  efficiency  of  the  eddy  transport 
processes  occurring  in  major  and  "minor"  breakdowns  of  the  polar  night  vortex. 
Here,  a  basic  assumption  is  that  the  ability  of  a  given  circulation  event  to 
produce  an  effective  (say  northward)  transport  of  trace  substance  can  be 
indirectly  inferred  without  any  detailed  knowledge  of  the  distribution  of  the 
substance  itself.  This,  of  course,  is  a  questionable  assumption,  but  has 
been  a  necessary  one  in  the  absence  of  any  knowledge  about  the  details  of  the 
trace  substance  distribution. 

Article  D  is  the  result  of  a  diagnostic  study  which  attempts  to  determine 
the  mean  transverse  circulation  about  the  polar  front  jet  stream.  Previous 
studies  have  shown  that  significant  amounts  of  stratospheric  air  can  enter  the 
troposphere  through  the  polar  front  Jet  stream  "front"  in  association  with 
specific  storm  systems.  Such  studies,  however,  do  not  answer  the  question 
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of  whether  there  is  a  weak  (but  systematic)  exchange  of  air  between  the 
stratosphere  and  troposphere,  resulting  from  a  mean  transverse  circulation 
about  this  jet  stream  system. 

In  this  study  such  a  calculation  was  performed,  using  vertical  motions 
from  a  "complete"  (Rossby  number  <  1)  diagnostic  balance  adequation  model. 
These  derived  vertical  motions  were  averaged  with  respect  to  the  polar  front 
jet  stream  core.  The  results  indicated  that  a  direct  mean  transverse  circu¬ 
lation  exists  about  the  jet  core  and  may  indicate  a  slow  import  of  stratospheric 
air  into  the  troposphere.  This  transverse  circulation  also  was  shown  to  give 
consistent  results  in  the  mean  kinetic  energy  balance.  Thus,  such  a  diagnostic 
balance ,  adequation  model  appears  to  be  a  very  useful  tool  for  conducting 
the  so-called  "dynamical  process"  studies. 

Articles  E  and  F  are  also  "traditional "'studies  in  that  they  attempt  to  infer 
transport  processes  without  detailed  knowledge  of  the  distribution  of  the 
trace  substance  (in  these  cases,  radioactivity)  in  the  atmosphere.  However, 
here,  as  in  a  number  of  previous  studies  by  other  authors,  information  is 
available  on  the  distribution  of  radioactivity  at  the  earth's  surface.  This  is 
very  valuable  because  the  "history"  of  a  given  trace  substance  can  often  be 
traced  from  source  to  sink  when  supplemented  by  detailed  meteorological 
analysis.  Article  E  falls  into  this  category. 

The  usefulness  of  accompanying  surface  radioactivity  is  shown  in  Article 
F.  From  previous  individual  case  studies,  import  of  stratospheric  air  into 
the  troposphere  in  the  polar  front  jet  stream  region  has  been  related  to  strong 
cyclogenesis  is  a  circulation  feature  which  should  serve  as  an  indication  of 
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startosphenc-tropospheric  mass  exchange5.  Article  T  is  an  effort  to  relate 
an  inoz-x  which  is  an  indicator  of  cyciogeriesis  to  the  radioactivity  fluctu¬ 
ations  at  the  ground.  The  results  show  that:  there  exists  a  significant 
cc.relatiori  between  the  two — thus,  establishing  the  mechanism  on  a  firm 
basis.  Without  the  surface  radioactivity  measurements,  however,  such  an 
inference  would  have  been  difficult,  if  not  impossible. 

Articles  G  and  H  both  employ  "new"  techniques  to  investigate  trace  sub¬ 
stance  transport  problems.  They  both  avoid  or  reduce  the  problem  of  inade¬ 
quate  trace  substance  data  resolution  by  employing  specific  techniques. 

Article  G  is  an  attempt  to  provide  direct  quantitative  calculations  of  strato¬ 
spheric  ozone  transport,  using  data  from  the  Air  Force  Cambridge  Research 
Laboratories  North  American  Ozonesonde  network.  In  the  past  this  network 
has  been  very  difficult  to  utilize  for  large-scale  transport  studies  even 
though  the  individual  observations  have  proven  to  he  quite  reliable.  The 
basic  problem  with  the  data  has  simply  been  one  of  inadequate  resolution. 
Observations  are  normally  taken  only  once  weekly  in  a  network  of  just  eleven 
stations  over  the  North  American  continent.  Such  a  network  is  too  sparse  to 
be  able  to  routinely  establish  synoptic  analyses  of  ozone  mixing  ratio  which 
possess  spatial  and  temporal  continuity. 

In  this  study  the  time  resolution  problem  was  minimized  by  considering 
only  a  specific  period  in  early  May  1963  in  which  the  ozonesonde  observations 
were  taken  daily.  Also,  the  effective  horizontal  data  resolution  was  increased 
by  computing  24-hour  forwaid  and  backward  isentropic  trajectories  from  each 
of  the  ozonesonde  observation  points  and  assuming  conservation  of  ozone 
mixing  ratio.  Ideally,  such  a  technique  increases  the  effective  number  of 
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data  points  by  a  factor  of  three.  In  practice,  because  the  trajectories 
leave  the  network  region  and  all  stations  did  not  report  every  day,  the 
number  is  somewhat  less.  Use  of  this  technique  does  allow  ozone-mixing 
ratio  analyses  to  be  prepared  which  are  consistent  in  time  and  space. 

With  such  analyses,  it  then  becomes  possible  to  compute  ozone  trans¬ 
port  quantities  directly.  In  this  study,  the  calculations  were  performed  by 
integrating  the  ozone  continuity  equation  over  an  arbitrary  volume  and 
separating  the  flux  convergence  terms  into  "mean"  and  "eddy"  effects, 
relative  to  the  computational  volume.  This  calculation  was  quite  successful 
in  that  it  showed  the  existence  of  a  systematic  downward  flux  of  ozone  due 
to  synoptic-scale  vertical  motions .  The  vertical  flux  by  the  area-mean 
vertical  motions  was  found  to  be  quite  sensitive  to  small  errors  in  the 
specification  of  the  diabatic  heating.  In  future  studies  of  this  kind, 
considerable  care  may  have  to  be  exercies  in  order  to  reduce  this  difficulty. 

In  Article  H,  an  attempt  is  made  to  simulate  numerically  the  transfer  of 
trace  substances  in  the  stratosphere.  In  this  type  of  approach,  the  problem 
of  lacking  trace  substance  data  is  circumvented  because  the  distribution 
of  the  trace  substance  is  evolved  numerically  from  an  arbitrary  initial  state. 

Previously,  such  numerical  approaches  have  been  separated  into  two 
distinct  categories.  The  first  and  simplest  method  is  to  treat  all  the  possible 
transport  mechanisms  as  a  diffusion  process.  The  second  approach  is  to 
introduce  a  tracer  into  a  numerical  model  of  the  general  circulation  and  study 
the  subsequent  evolution  of  the  tracers.  This  is  potentially  a  very  powerful 
method  which  will  see  considerable  application  in  the  future.  The  first  method 
is  limited  because  of  the  very  arbitrary  specification  of  the  diffusion  process. 
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The  second  method  is  only  limited  in  that  a  given  general  circulation  model 
may  or  may  not  properly  simulate  the  real  atmosphere  with  sufficient  fidelity 
to  reproduce  the  observed  transport  characteristics.  At  present  this  may  be 
particularly  true  in  the  region  of  the  stratosphere,  especially  in  regard  to 
seasonal  variations  in  the  circulation. 

Because  of  these  difficulties,  an  intermediate  approach  is  attempted  here. 
The  observed  winds  in  the  stratosphere  on  a  daily  basis  are  used  to  advect 
the  arbitrarily  chosen  tracers.  From  these  winds  and  the  temperature  fields, 
vertical  motions  are  obtained  on  the  same  scale  as  the  horizontal  wind  data. 
These  horizontal  and  vertical  winds  are  then  adjusted  slightly  so  that  they 
yield  a  consistent  mean  meridional  circulation  when  zonally  averaged. 

Therefore,  these  winds  give  a  quite  reliable  specification  of  the  atmospheric 
structure  for  motions  on  the  scale  of  the  meteorological  network  or  greater. 

The  motions  below  the  resolution  of  the  grid  are  treated  as  a  sub-grid-scale 
diffusion  process. 

The  Eulerian  form  of  the  continuity  equation  for  a  given  trace  substance 
is  then  integrated  numerically  in  time  on  a  three-dimensional  grid,  using 
the  derived  horizontal  and  vertical  motions  to  provide  the  necessary  advections. 
This  technique  gives  very  encouraging  results  in  that  the  evolved  tracer 
fields  very  shortly  begin  to  possess  characteristics  quite  similar  to  the 
limited  knowledge,  available  about  the  observed  fields. 

This  result  suggests  that  such  a  technique  might  be  applied  to  a  number 
of  specific  transport  problems  for  different  types  and  scales  of  circulations  in 
both  the  troposphere  and  the  stratosphere. 
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This  is  a  compilation  of  research  results  performed  under  AEC  Research 
Contract  AT(49-7)-3206.  All  of  the  various  efforts  have  been  directed 
toward  gaining  a  better  understanding  of  the  atmospheric  processes  which 
lead  to  systematic  large-scale  trace  substance  transport.  A  variety  of 
approaches  have  been  employed  in  this  and- in  earlier  work  on  this  contract. 
The  last  article  in  this  report  provides  an  overview  and  critique  of  the 
applicability  of  these  various  methods. 


